﻿J4 , Lol UNIVERSITATEA Ш CLUI-HAPOCA Facultatea de Matematică si Fizici и T OPROIU Â PÂl V POP V URECHE CULEGERE DE EXERCIȚII, PROBLEME Șl PROGRAME DE CALCUL I 1909 UNIVERSITATEA DIN CLUJ-NAI-OCA FACULTATEA DE МАТШАТІСА ȘI FIZICĂ T Or ROI U Ă PAL V, POP V URECHE ASTRONOMIE CULEGERE DE EXERCIȚII, PROBLEME SI PROGRAME DE 9 CALCUL EDIȚIA А II-A f CLUJ-КАРОСA 1989 - 2 - гаоиНАМШ 4: Determinare® datei iuliane ceraspunzătaare anei datо calandariatiea (din oalaadarul gregorian ■au iulian) • 262 PROQRAJtpii’1 г Transformarea coordonatelor orizontale îri ■ coordonate', ecuatoriale', ;invers-' « ț265 PROGRAMUL 2: Transformarea coordonatelor ecuatoriale ■ date' pentru spoea 1950 0 în coordonate ecuatori- ; pentru epoca t s ținîndu-se seama de prece- sie 269 PROGRAMUL Зг Calcularea coordonatelor geocentrice Ecuatoriale rectangulare, ale unui punct de pe suprafața Pămîniului pentru care se cunosc coordonatele geodezice 272 - 4 - precum și я unor tabele, e căror: analiză conduce la descifrarea unor legități în studiul fenomenelor astronomice Spre deosebire de lucrări cu caracter riguros de cercetare științifică, in lucrarea de față precizia calculelor numerice nu este esențială, La aceea, am admis obținerea soluțiilor cu diferite precizii și sîntam de părere că și cei oe vor folosi prezenta culegere vor putea rotunji datele exacte din exercițiile și problemele formulate, firește, în interiorul unor limite raționale,' Această ultimă observație-nu se referă însă la programele de calcul, a căro-y includers în prezenta culegere are și scopul de' a venit'îrv'ajutorul celor-care -se inițiază în cercetarea astronomică, în cadrul cercului științific studențesc sau al elaborării lucrărilor de diplomă Programele sînt elaborate pentru microcalculatoare, în limbajul BASIC Rezultat al unei colaborări,colegiale de mai mulți ani, volumul a fost redactat astfel; coordonarea și capitolele III, VI, VII - de prof dr A Păi , capitolele VIII-XV - de lect- dr V Ureche , capitolele J, II, IV, V - de lect dr V Pop, iar programele da calcul - de cerce t șt dr T Oproiu ■Autorii mulțumesc, și pe această cale cercetătorului principal dr І Todoran, cercetătorului științific dr V ¥ioc și analistului programator B Pârv pentru valoroasele-lor observații și sugestii exprimate în vederea îmbunătățirii conținu-tului lucrării ' ; AUTORII - 3 - ? Й S ? А T I Exercițiile și problemele aînt menite să faciliteze însușirea cursului predat, reprezentând aprofundări ale bazelor teoretice ale astronomiei, iar programele de calcul înlesnesc efectuarea unor calcule mai laborioase specifice acestei discipline Lucrarea de față urmează ordinea și conținutul capitolelor manualului publicat de A Păi și V Ureche: Astronomie (Editura didactică și pedagogică, București, 1983), adresat studenților facultăților sau secțiilor de matematică ale u-•niversitățilcr din țară, potrivit cu planul de învățămînt în vigoare Excepție face capitolul I, care cuprinde noțiuni și formule necesare în unele capitole, din geometria pe sfera și trigonometrie sferică Plecare capitol al culegerii conține exerciții și probleme rezolvate de autori, pentru a explica, pfin exemplele alese, conținutul mai dificil al chestiunilor teoretice din manual; este, de asemenea, propus pentru rezolvare individuală un număr minim de probleme, la care am dat numai răspunsu- iar în unele sasuri și indicații pentru soluționarea lor Volumul cunoștințelor înglobate în exemplele noastre diferă de la un capitol la altul, aceasta și în funcție de gradul de completitudine a explicațiilor conținute în diversele părți ale manualului la care tie referim Autorii au socotit utilă prezentarea unor figuri matematice, care își au rolul cunoscut în rezolvarea problemelor, PARTEA î ВТ î I EXERCIȚII ȘI RROBLEME - 6 - CAPITCLU ELEMENTE DE GEOMETRIE PE SPERI ȘI DE TRIGONOMETRIE SFERICĂ 1 1 Care sînt elementele de bază ale geometriei sferice? Răspuns : Geometria sferică se ocupă cu studiul figurile, de pe o sferă Dreptelor și segmentelor de dreaptă din geome- tria plană li rrespund cercuri și arce de cercuri pe sferă Deosebim' cercuri жагі, ale căror plane trec prin centrul sferei; - cercuri mici căror plane nu trec prin centrul sfe- rei Prin două puncte ale sferei trece un cerc mare și numai unul, dacă cele două puncte nu aînt diametral opuse în ade- văr, prin două puncte А, В de pe sferă și prin centrul O tre- ce un plsn unic dacă punctele A, 0 și В nu eînt coliniare In- tersecția acestui plan cu sfera determină cercul mare ce tre- ce prin A și В (fig 1 1) Fig 1 1 Cercuri pe sferă Diametrul sferei per- pendicular pe planul unui cere intersectează suprafa- ța sferei în două puncte nu- mite polii cercului Două cercuri- mari de pe aceeași sferă se intersectează întotdeauna în do- uă puncte diametral opuse - 7 - Unghiul format de două arce de pe sferă, care se intersectează într-un punct, este, prin, definiție, unghiul (s, 180'') > o , de tangentele duse la cele două arce în punctul lor de intersecție In particular, măsura unghiului sferic APE (fig 1 1) este egală cu: - măsura unghiului diedru al planelor PAP' și PBP'; - măsura arcului AB de pe cercul mare ai cărui poli sînt P și P' Măsura unui arc de cerc mic de pe sferă se poate exprima cu ajutorul măsurii arcului de cerc mare Pentru a arata a-cest lucru, să considerăm cercul mic cu centrul în C și care are planul paralel cu planul cercului AB cu centrul în 0 Cer curile mari PAP' și PBP' determină pe cercul mic arcul DE Avem: măs(DE)=CD-măs(DCE) ; măs(AB) =OA-măs(ÂOB) ==> măs(DE)/măs(AB) = — = = cos S | OA măs(DCE) = măs(AOB) j și atunci putem scrie: măs(DE) = măs(AB) • cos S Fig 1 2 Fusul sferic Fusul sferic este figura de pe sferă formată din două semicercuri mari ale căror extremități coincid cu extremitățile diametrului lor comun (fig •1 2) Elementele fusului sferic sînt: - două laturi, ABA' și АСА', de 180° fiecare; в - două unghiuri congruenta, A și A' Dată fiind măsura unghiului A al fusului sferic de pe sfera de rază R, se poate determina aria fusului, utilizind regula de trei simplă, și anume: la unghiul cu măsura 25Г со- 2 л respun de aria 4-TtR , iar la măe(A) corespunde aria rezultă: SA = (4%R2/(2Jt)) • măa(Â) = 2R2măs(Â) , (1 2) care exprimă faptul că aria fusului sferic este proporțională cu măsura unghiului său Poligonul sferic este figura de pe sferă mărginită de arce ie cerc mare (laturile poligonului) mai mici de ISO0 și limitate de intersecțiile lor consecutive Un poligon sferic se numește convex dacă este situat de aceeași parte în raport cu fiecare din cercurile mari cărora le aparțin laturile sale £n caz contrar, poligonul se numește concav 1 2t Ce este triunghiul sferic? Răspuns ; Triunghiul sferic este cel mai simplu poligon sferic convex Elementele triunghiului sferic sînt: trei laturi a, b,c și trei unghiuri А, В, C {fig 1 3), fiecare din acestea fiind mai mici de 180c (triunghiuri euleriene) Dacă un unghi al triunghiului sferic este de 90°, triunghiul se numește dreptunghic Triunghiul sferic care are o latură de 90° se numește rectilater Există triunghiuri sferice tridreptunghice și trirectilatere (exemplu: triunghiul format de ecuator și două meridiane perpendiculare între ele) - 9 - Fig 1 3 Triunghiul sferic Fig 1 4 Triunghiul polar al unui triunghi sferic dat Unind vîrftJ'ile л, В, C ale triunghiului sferic ABC cu centrul 0 al sferei de rază R, se obține triedrul OABC corespunzător triunghiului sferic ABC între fețele triedrului și laturile triunghiului sferic există relațiile: măe(AOB) = măs(c) , măa(BOC) = măs(â) , măs(ACC) = măs(b) Un unghi al triunghiului sferic, de exemplu cel din A, este definit de tangentele duse în A la arcele de cerc AB și AC Deoarece aceste tangente sînt perpendiculare pe OA, rezultă că ele formează tocmai unghiul plan al diedrului definit de feței-ДОС și AOB In baza proprietăților generale ale triedrului, deducem că: ■ ' 1) O latură a triunghiului sferic este mai mică decît sumr; celorlalte două; diferența a două laturi este mai mică decît latura a treia: IC в «> / B’ « 180°-b (1 6) c = 180°- C'J [ C’ = ISO0- o Din relațiile (1 5) și (1 6) rezultă că dacă triunghiul dat este dreptunghic, atunci triunghiul său polar este recti-later și invers Pornind de la inegalitatea (1 3) scrisă pentru triunghiul sferic A'B'C’: a'4b' + o* și ținînd seama de relațiile de tipul (1 5), se obține 180°- A 0°), a-dică suma unghiurilor triunghiului sferic ABC este mai mare de-cît 180° ! 1 4 Să se deducă formula pentru aria triunghiului sferic ABC i , ; - Rezol v a r e : Aria triunghiului sferic ABC de pe sfera de rază R (fig 1 5) se poate c ded uce ușor, dacă observăm că surea -Ariilor celor trei fuse/Sferice (AA1 , BB' , CC') care eu drept unghiuri un- ghiurlle triunghiului sfe- ric dat acoperă emisfera vizibilă a sferei piue de doua ori aria a triunghiului sferic dat Ținînd Plg 1 5 Pentru determinarea ariei triunghiului sferic seama și de faptul că aria - 12 fusului sferic este dată de (1 2), rezultă: 2R2măs(Â)+2R2măe(B) + 2R2măs(C) - 2Sabc = 2KR2 , О Л Л Л sau: R (măs(A)+măs(B)+măs(C) - 1Г) “ $ABC de unde, utili- zînd notația (1 9), se obține: SABC = Er2 sau SABC = №r2/180°)E • (1 10) Dacă R = l, rezultă că E = sabq • 1 5 Să se deducă formulele ce exprimă teoremele cosinusului, sinusului și formula celor cinci elemente referitoare la laturile triunghiului sferic ABC (formulele lui Gauss) Rezolvare : Se consideră triunghiul Mg 1 6 Pentru deducerea formulelor de baza ale trigonometrie! sferice sferic ABC de pe sfera cu centrul în 0 și de rază unitate (R=l) Alegem sistemul trirectangu-lar de referință Oxyz astfel încît axa Oz să intersecteze sfera în’ vîr-ful A al triunghiului, iar vîrful В să aparțină planului xOz (fig 1 6) Componentele vectorilor de poziție ai vîrfu- - 13 - rilor triunghiului eînt: (0, 0, 1) ; '?5 (sine, 0, cose) ; e*- (sin b cos A, sinbsinA, coeb) Făcînd produsul scalar al vectorilor unitari e*,- si el obținem: (ё^, е\) - В е^й’И e^fi cosCe’g,е^) = cos b сов c + sin b sin c cos A , sau: cos a « cos b cos c + sin b sin c cos A (1 11) V ■ -•■■■ у în mod analog se obțin și formulele: cos b « cos а сов o + sin a sin c cos В , (1 11 1) cose = cos a cos b + sin a sin b сов C (1 11") Aceste formule reprezintă teorema cosinusului referitoare la laturile triunghiului sferic ABC înmulțind (1 11) cu cose și adunînd membru cu membru la relația (1 11')> efectuînd reducerile și împărțind prin sine relația rezultată, se obține: sin a cos В x cos b sin c - sin b сов c cos A ) (1 12) în mod analog se obțin încă două relații ce reprezintă, împreună cu (1 12), formulele celor cinci elemente referitoare la laturile unui triunghi sferic ABC Din (1 11) rezultă: cos A = (cos a - coe b cos c)/(sin b sin c) Ridicînd această relație la pătrat, scăzînd fiecare membru al m 2 egalității din 1 și împărțind prin sin a fiecare membru al relației rezultate, se obține': 2 2 2*2 2 sin A/sin a x (1 - cos a - cos^b - сое c + + 2 cos a cos b сов c J/(sin‘ а віп*"Ъ sin£c) Avînd în vedere faptul că membrul din dreapta al acestei reia- 14 ții este o funcție simetrică în raport cu a, b, c, rezultă că; sin A/sin a = sin B/sin b sin C/sin c (1 13) Aceste formule reprezintă teorema sinusului referitoare la laturile unui triunghi sferic ABC Relațiile: cos a = cos b сое c + sin b sin с сое A sinacosB = cos b sin c - sin b cos c cos A , (1 14) sin a sin В = sin b sin A constituie formulele fundamentale ale trigonometrie! sferice sl poartă numele de formulele lui Gauss Prin permutarea circulară a literelor din (1 14) se obțin trei grupe de relații fundamentale: cos a = cos b cos c + sin b sin c cos к , (I) cos b = cos c cos a + sin c sin а сов В cos c = cos a cos b + sin a sin b cos 0 ; sin a cos В = cos b sin c - sin b рое c cos A , sin a cos C = cos c sin b - sin c cos b cos A , sin b cos C = cos c sin a - sin c cos a cos В , (II) sin b cos A = cos a sin c - sin a cos c cos В , sin c cos A = cos a sin b - sin a cos b cos C , ein c cos В = сое b sin a - sin b сое а сов C ; sin a sin В = sin b sin A, (III) sin b sin C = sin c sin В , sin c sin A = sin a sin C Prin Introducerea elementelor auxiliare m și M, date de relațiile- - 15 - tgК * tg b cos A , (1 15) m = сов b/сов И = sin b cos A/ein M (deduse din notațiile: сов b = m сов И și sin Ъ сов A = m ein M), primele două formule din (1 14) se transformă în două expresii, calculabile prin logaritmi, și anume: cos s = m cos(c - M) , (1 16) sin a cos В = msin(c-M) 1 6 Să se deducă formulele lui Gauss referitoare la un-ghiurile triunghiului sferic ABC Rezolvare : Scriind formulele lui Gauss (1 14) pentru triunghiul polar A’B'C' al triunghiului ABC și ținînd seama de relațiile (1 5) și (1 6), se obține: сое A = - cos В cos C + sin В віп С сов а , sin А cos b = cos В sin С + сов С sin В сов а , (1 17) \ ' ' ■ • ? g sin А sin b sin a ein В , adică formulele lui Gauss referitoare la unghiurile triunghiului sferic ABC Observați e : Formulele lui Gauss pot fi deduse și prin, alte metode ca: metoda rotației axelor de coordonate, metoda matricială, metoda vectorială etc Să se particularizeze formulele lui Gauss pentru triunghiul sferic dreptunghic (A = 90°) - 16 - Răspuns: сова = cos b cos с = ctgВctgC sin b =• віп a sin В = tg c ctg C , sin c = sin a sin C = tg b ctg В , cos В = cos b sin C = tg c ctg a , cos C '=' cos c sin В = tg b ctg a 1 8 Să se particularizeze formulele lui Gauss pentru tri-~ unghiul sferic rectilater (a = 90°) Răspuns : ' cos A - сов В сов = - ctg b ctg c sin В sin A sin b tg C ctg c sin C sin A sin c tg В ctg b cos b сое В sin c - tg C ctg A cos c cos C sin b - tg В ctg A C 1 9 j Să se deducă formulele pentru cos(A/2) , sin(A/2) , ~ cos((A-B)/2), sin((A-B)/2), cos((A+B)/2) sin((A+B)/2) ! în funcție de laturi, folosind notația a+b + c = 2p (for-I mulele lui Delambre) Re salvare : Din formula cosinusului deducem ușor: cos A = (cos a - cos b cos c J/țsin b sin c ) Răcind combinațiile 1 + cos A = 2 cos (A/2) și analog 1 - cos A = 2 = 2 sin (A/2), obținem pe A/2 cu ajutorul formulei; 1 + cos A = (cos a - cos(b+c) )/(sin b sin c ) = = 2 sin p sin(p-a)/(sin b sin c) , unde am introdus notația a+b+c=2p - 17 - Atunci obținem imediat: coe(A/2)= + Veiiip sin(p-a)/(sin b ein c) , (1 18) deoarece A 0 în mod șnalog ѳе obține formula: sin (A/2) = + \fsin(p-b) sin(p-c)/(einb ein c) (1 19) și, prin permutări circulare, formule analoge pentru unghiurile B/2 și C/2 Din formulele (1 18), (1 19) și din formulele similare pentru celelalte unghiuri ее deduce ușor: coe(A/2) cos(B/2) = (sin p/sin c) sin(C/2) , sin(A/2) sin(B/2) = (sin(p-c)/Bin c) ein(C/2) , de unde, adunînd membru cu membru și ținînd seama de faptul că 2p = a + b +■ c , obținem: cos((A~B)/2) = (sin((a+b)/2)/sin(c/2)) sin(C/2) (1 20) Formulele, pentru cos ((A+B)/2), sin((A-B)/2) și sin((A+B)/2) se obțin în mod analog Acestea reprezintă formulele lui Delambre, care se scriu de regulă sub forma: cos((A-B)/2)/sin(C/2) = sin((a+b)/2)/sin(c/2) , cos((A+B)/2)/sin(C/2) = cos((a+b)/2)/cos(c/2) , (1 21) sin((A-B)/2)/cos(C/2) = sin((a-b)/2)/sin(c/2) , sin((A+B)/2)/coe(C/2) ® cos((a-b)/2)/coe(c/2) 1 10 Să se deducă formulele pentru tg((A-B)/2)/ctg(C/2) tg((A+B)/2)/tg(C/2) în funcție de laturi și formulele pentru tg((a-b)/2)/tg(c/2) , tg((a+b )/2 )/tg(c/2) în ! 18 ' funcție de unghiuri (formulele lui Reper) i — i -1 Re zo Ітат e : Din formulele lui Delambre (1 21), prin împărțire, obținem: tg((A-B)/2)/ctg(C/2) = sin((a-bj/2)/sin((a+b)/2) , tg((A+B)/2)/tg(C/2) = coe((a-b)/2)/cos((a+b)/2) , (1 22) tg((a-b)/2)/tg(c/2) = віп((А-Б)/2)/віп((А+В)/2) , tg((a+b)/2)/tg(c/2) = cos((A-B)/2)/cce((A+B)/2) , care HÎnt formulele (sau analogiile) lui Neper (sau Eapier) ; 1 11 Să se exprime excesul sferic £ în funcție de laturile triunghiului sferic (formula lui D’Huillier) Rezolvare : Vom considera și utiliza următoarea formulă din trigonometrie plană: tg((x-y)/2) = (sinx - siny)/(cos x + cos y) Deoarece avem: tg(e/4)«tg((A+B+C-%)/4) = tg(((A+B)/2 - (1T-C)/2)/2) , pentru x = (A+B)/2 și у = (7Г-С)/2 , rezultă ușor: tg(£/4) = (віп((А+В)/2) - сов(С/2))/(сов((А+В)/2) +sin(C/2)) înlocuind în această relație sin((A+B)/2) și cos ((A+B )/2) prin expresiile date de formulele lui Delambre (1 21), obținem după calcule simple: tg(£/4) = (2 sin((a-b+c )/4) sin((c-a+b)/4 )/(2 coa((a+b+ +c)/4) cos((a+b-c)/4)))(l/tg(C/2)) Exprimînd tg(C/2) cu ajutorul formulei corespunzătoare ce se obține din (1 18) și (1 19), relația de mai sus devine: - 19 - tg(E/4) = (ein((p-a;/2) ein((p-b)/2)/(cos(p/2) cos((p-c)/2))) • • Vsîn p 8Іп(р-с)/(віп(р-а) sin(p-b)) Ridicînd le pătrat și folosind formula din trigonometrie plană sini - 2 ein(x/2) cos(x/2), rezultă: tg2(£/4) = tg(p/2) tg((p-a)/2) tg((p-b)/2) tg((p-c)/2) , (1 23) care este formula lui L'Huillier 1 12 I Să se arate că în cazul limită, cînd R—» oo și tri- i unghiul sferic se apropie de un triunghi plan, teoremele sinusului și cosinusului ale trigonometrie! sferice ; ее transformă în teoremele corespunzătoare din trigonome- j 1 i tria plană, iar din formula lui L’Huillier rezultă formula lui Heron din trigonometrie plană Rezolvare : Să considerăm trei puncte din spațiu ce nu sînt coliniare, А,В, C Ele determină un plan și, în acest plan, un triunghi în același timp, există o infinitate de sfere pe a căror suprafață se găsesc cele trei puncte Dacă aceste sfere sînt ordonate după mărimea razelor, atunci cînd R-pgd curburile respective descresc și triunghiul sferic se apropie de un triunghi plan; unghiurile de pe sferă se aprc pie de unghiuri plane obișnuite și excesul sferic E devin, oricît de mic Lungimile laturilor triunghiului plan, a, b au drept corespondente în triunghiul sferic arcele a/R, b/R, c/R, măsurate în unități de arc Din dezvoltarea funcțiilor trigonometrice sinus și cosinus în serii convergente, notînd ă/R = qa, b/R=qb, c/R=qc și prin Sj - mărimile ce tind la 0 odată cu 1/R^, se obține: 20 - s±nqa= Ча-фЗ! + -qa(l-qa/6+ 8J , 2 ' 1 cos qa = 1 - qfi/2 + S2 • Din teorema sinusurilor din trigonometrie sferică rezultă: sin oc : sin p : sin Y = ă(l - q2/6 + 8j) : b(l - q^/6 + 5-j) : : " qc/6 + S55 și la limită se obține: sin OC : sin p : sin Y = ă ,s b : c , (1 24) adică teorema sinusurilor din trigonometrie plană Din teorema cosinusurilor din trigonometria sferică se obține: 1 - q2/2 + S2 = (1 - q2/2 + 8Д)(1 - q2/2 + S6) + + " 4t/6 + S3)(l - q|/6 + Sș) cos oc , de unde: “ qa/2 + S2 = " (qb + %)/2 - ^4,6 + qb%(1 ~ o 2 - (q^ +Qc)/6 + £- 5) cos oc și la limită se obține: ~“9 —2 - — — a = b + c - 2 b c cos ot , (1 25) adică teorema cosinusurilor din trigonometria plană Dacă în formula lui L'Huillier (1 23) funcțiile tangentă se înlocuiesc cu arcele unghiurilor, acestea'din urmă fiind foarte mici, obținem formula aproximativă: e/4 = (1/4) \Г(?/Н) ((p-E)/R) ((p-b )/R) ((F-c )/K) ■ Utilizînd expresia de mai sus și ținînd seama de expresia (1 10) pentru aria unui triunghi sferic, se obține formula lui Hehcn pentru, triunghiul plan: 21 SĂBC ж £r2 = = Vp(p-a)(p-b)(p-c) (1 26) [ 1 13 |Să se arate că într-un triunghi sferic echilateral există relația: cosa = cos A/(l - cos A) Rezolvare : Din ipoteză rezultă a «= b = c în aceste condiții, prima formula a lui Gauss (1 14) devine: cose = coe^a + sinca cos A , din care deducem succesiv: O сое a (1 - cos a) = (l-cos a) cos A , cos a = (1 + cos a) cos A , сое a (1 - cos A) = cos A De aici rezultă relația căutată j 1 14 Să se determine aria triunghiului sferic echilateral Г” care are unghiurile de 75°• Rezolvare : Din ipoteză avem A = В = C = 75° = 53Г/12 Excesul sferic ‘este g = A + В + С -5Г = 3 • 51Г/12 - ЗГ - 7Г/4 Apli cînd formula (1 10), găsim aria triunghiului: SABC = E R2 = 7Îr2/4 • 1 15 Știind că aria unui triunghi' sferic dreptunghic isoscel este o treime din aria unui cerc mare al sferei, să se calculeze unghiurile triunghiului Rezolvare : Conform ipotezei, avem A = 90°, B = C Aria 22 - triunghiului sferic este atunci: S,™ « R2(0l/2 +B + B-ЭТ) = й2(2В-1Г/2) 2 Dar din ipoteză rezultă că această arie este egală cu 3tR /3 2 2 • Atunci putem acrie că R (2В-7Г/2) = 3TR /3, de unde rezulta; В = (7Г/3 + 3C/2)/2 = 53Г/12 , sau, în grade, В - С - 75° 1 16 ; Să se determine lungimea laturii unui triunghi sfe- — așezat pe o sferă cu raza R = 5m, știind că a = = 32°17’15” Rezolvare : Calculăm a în secunde de arc (a"): 32° = 32 • 3600" = 115 200" 17’ = 17 • 60" = 1020" 15" = 15" a" = 116 235" Dacă R este raza unui cerc și oc numărul de radiani al unui arc de pe acel cerc, atunci lungimea arcului este: l = Ro, (1 27) Să calculăm numărul de radiani al unui arc de 1", număr pe care-1 notam cu C0 Știind că unui arc de 180° (=■ 648 000") îi corespund 31 radiani, prin regula de trei simplă găsim: щ = Я /648 000 = 3,14159/648 000 = 1/206 265 rad Atunci, din (1 27), unde ot = a"co , calculăm imediat Iun- gimea laturii & : L = Ra"as = 5 - 116 235 • (1/206 265) ’= 2,817 m - 23 CAPITOLUL II ASTRONOMIA SFERICA 2 1 Enumerați punctele fundamentale de referință de pe sfera cerească Răspuns : Polul nord al lumii (P), polul sud al lumii (?’); zenitul (Z), nadirul (Z* ); punctul cardinal sud (S), punctul cardinal nord (H), punctul cardinal est (E), punctul cardinal vest (V); polul nord ecliptic (TT) , polul sud ecliptic (TT' ) ; punctul vernal (ф = punctul echinocțiului de primăvară), punctul autumnal Gfi= punctul echinocțiului de toamnă), punctul solstițiului de vară (*£ = Q = 90° = 6h , - S 1 m (pentru culminația superioară; zg, - distanța zenitală a astrulai; semnul minus se ia atunci cînd astral se află la sud de zenit, iar semnul plus - la nord de zenit); zm = ieo°- (ч> + s) (pentru culminația inferioară) Indicații ; Se consideră fig 2 5, unde am presupus că aștrii ee afla în spațiu la distanțe foarte mari de P&nînt, în-cît ei pot fi considerați ficși pentru timpul perioadei de rotație a Pămîntului Am utilizat următoarele notații: pp' - axa de rotație a Pămîntului (T - centrul Pământului)-qq' - ecuatorul Pămîntului, săgeata indicînd sensul rotației terestre; 0 - poziție observatorului pe meridianul terestru pqp' ; OZ - linia verticalei locului 0; KS - meridiane locului 0 (intersecția orizontului matematic cu planul figurii); OP - axa lumii; 28 QQ'f - intersecția planului ecuatorului ceresc cu planul figurii (OPlIpp', QQ" || qq'); ~e vede că planul figurii eate planul meridianului ceresc și că planul meridianului ceresc coincide cu planul meridianului geografic; O-, - poziția observatorului după o jumătate de perioadă de rotație a Pămîntului; - linia verticalei locului' 0- ; NfSi “ poziția orizontului locului Сц : 0-, P - axa lumii; Planul ecuatorului ceresc QQ* răniîne neschimbat H g 2 5 Pentru stabilirea condițiilor de vizibilitate a aștrilor Planul orizontului observatorului descrie o suprafață conică, unghiul de la vîrf al conului de sus fiind , iar - 29 - unghiul de la vîrf al conului de jos - SKS-, Din fig 2 5 se vede că toți aștrii care ee află în conul de sus sînt întotdeauna deasupra orizontului (deci nu apun), iar aștrii care se află în conul de jos se află întotdeauna sub orizont (deci nu răsar niciodată pentru observatorul dat) Pentru deducerea formulelor pentru distanța zenitală în cazul trecerii la meridianul ceresc al locului, am indicat în fig 2 5 direcțiile către trei stele: Vega, Sirius și % Dre , care se află în spațiu aproximativ în același plan (în fig 2 5, exact în planul figurii) F2 șj ' Să se descrie fenomenele~mișcării aparente a Soarelui j pe sfera cerească, pentru următoarele latitudini geografice: 4>=0°; +2J°,5î +66°,5; +90° (din emisfera boreală) Cum ee desfășoară fenomenele similare pentru emisfera aus-relă ? Rezolvare : Vom examina pe rînd cele patru cazuri menționate în enunț, în proiecția sferei cerești pe planul meridianului ceresc al locului de observație (0) I Cazul^j^O^lgcuațorul,} (fig 2 6) Axa lumii (PP*) coincide cu meridiane (NS): polul nord al lumii (P) coincide cu punctul cardinal nord (N), iar polul nud al lumii (P') coincide cU punctul cardinal sud (S); ecuatorul ceresc (QQ') este perpendicular pe orizont (KS), toți paralelii cerești sînt perpendiculari;pe orizont și sînt îm-părțiți în cîte două părți egale; drept urmare, drumurile tuturor aștrilor deasupra orizontului, respectiv dedesubtul o-rizontului, sînt egale în cazul mișcării Soarelui: ziua este - 30 - întotdeauna egală cu noaptea, Soarele aflîndu-se de două ori, în zilele echinocțiilpr, la miezul zilei, la zenitul observatorului (Z) Se poate întîmpla, însă, ca la 21 martie și 23 septembrie centrul discului solar să nu ее afle exact în zenit Pig 2 6 Pentru determinarea mișcării aparente a Soarelui coreapunzînd unui loc situat pe ecuatorul terestru Pig 2 7 Pentru determinarea mișcării aparente a Soarelui corespunzînd unui loc situat pe unul din tropice Fig 2 9 Pentru determinarea mișcării aparente a Soarelui corespunzînd unui loc eituat într-unul din polii tereștri Pig 2 8 Pentru determinarea mișcării aparente a Soarelui corespunzînd unui loc situat pe unul din cercurile polare - 31 II Саги1 ^₽= 23j 5 H £țropicul Racului) (fig 2 7) Axa lumii este înclinată pe orizont cu 23° 5 și, drept urmare, ecuatorul ceresc (QQ’) este înclinat pe linia verticală (2Z’) cu același unghi De aici rezultă că tropicul Racului trece prin zenitul observatorului (2), iar tropicul Capricornului - prin nadir (Z ’ ) Dacă este vorba de Soare, a-cesta se află pe tropicul Racului o singură dată pe an, în ziua solstițiului de varăj în această zi, la miezul zilei, Soarele se află la zenit (Z) Daca observatorul (0) se află la latitudini >23°,5 (în emisfera boreală), atunci axa lumii (PP!) se ridică deasupra orizontului, ecuatorul ceresc (QQ’), dimpotrivă, coboară spre orizont, iar tropicul Racului va trece la sud de zenit și, ca urmare, Soarele, în mișcarea sa, nu va mai ajunge la zenit Dimpotrivă, dacă *•? ) = В , rezultă: tg 2 sinW^ +M) / cos M , unde К se determină după formula (2 3) 2 B i Să se scrie formula pentru distanța dintre două puncte de pe sfera cerească, de coordonate ecuatoriale e'^Gx -L S^), 6^2 (c* 2’^2 )• Rezolvare : Să considerăm fig 2 15 Aplicînd în triun- ghiul sferic Rt), 6"7 formula cosinusului, obținem: x'= c- 62 “ arccos(sin S ■sin Sg + cos S-j cos 8g • •cos(oc2 - ot^)) Transformarea formulei Pig 2 15 Determinarea dis- tanței dintre două puncte de pe sfera cerească de mai sus într-o formă cal- culabilă prin logaritmi se poate efectua într-un-mod a- nalog celui prezentat în problema anterioară 2 9 Un vapor pleacă în ziua de 8 iulie, ora 6 dimineața, ’ " din portul B (fn = 48°23 ’ N, Lj = 4°3O' V), pe un cerc mare, spre portul 0 (4>2 = 5°17' N, L2 = 52°33' V) Ciad ajunge vaporul la destinație dacă viteza ea medie de deplasare este de 18 mile marine pe oră ? Rezolvare : Se consideră triunghiul sferic PBC (din fig 2 16), în care: PB = 90°- ’fj , PC = 90°- 4>2 , BPC = L2 - Ly Pentru a determina latura BC a acestui triunghi sferic, a- plicăm prima formulă a lui Gauss (formula/cosinusului): Big 2 16 Pentru determinarea distanței între punctele te- restre В și cos ВС = cos BP cos CP + + sinBP sin CP cos BPC , formulă care se mai poate scrie în modul următor: cos BC =■ sin Ч’і sin 4g + cos 4^ cos % cos (Lg - ) be aici, ținînd seama de datele problemei, se obține: BC = 59°16',6 , BC ceea, ce echivalează cu: BC = 3556,6 mile marine (1 milă marină = lungimea arcului de cerc mare de 1* de pe и ("era terestră) • Timpul necesar pentru ca vaporul să parcurgă această ii istanță este: - 42 - t- 3556,6/lB» 197^6 «В bile și 5h,6 = 8 sile și 6to Așadar, vaporul ye ajunge în portul O la date de 16 iulie , ora 12 2 10 j La ce distanță se găsesc orașele Paris = - =48°50'K, Lx-2°20’E), Berlin (= 52°31 ’ B, L- « 13°24' E) și București ( , PS = 90°-S , ZPS = H, II- = 180°-A Aplicînd formulele lui Gauss triunghiului sferic IZ- (în care a S 90°-S , b = z), se obține: cos (90°- 8 ) ~ cos z cos (90C- + віпЧ sin z cos A , cos 8 (sin e cos cc - cos 6 sin «,) = sin z sin A înmulțind a doua relație cu cos6,'iar pe a treia cu sin © și adunînd, obținem: cos S coa cb ~ cos z сов $ сое Ѳ + sin z sin (H, 8 ) —» (A,z) Prima transformare se realizează pe baza formulei H - © - ot (în cazul nostru, H - 1^44m36®5) Pentru ■ »r de-a doua transformare, se consideră trl«ng i-iul de pozi- 1 le PZS (fig 2 18), căruia i se aplică formulele lui Gauss (luînd as z , bsSO^-Sj Se obține: cos z = ein 8 sin Ф + cos X сов cos H - ЯІП z COS A a ein S сов Ф -sin z ein A = sinH сов ' , sin ‘-P cos S сов H (2 7) 4 nu ниЬ formă matricială ninzsin A -ein z cos A cos z 0 -ein Ч’ cos Ч5 ‘ 0 cos ein сов S cos 8 sin 8 sin H сов H (2 в) 1 0 0 4> Înlocuind în formulele (2 7) valorile date Ч\ H și 8 , se cb- - 46 - țin coordor- ele orizontale căutate: А = 55°43’2Ок z = 28°17'4O" 2 13 La o epocă dată steaua Capella are coordonatele ecuatoriale ( sin > = sin 8 sin £ + cos & cos 6 sin oc , (2 9) сое в cos Ti = cos oc cos 8 , - 47 - sau, sub formă matricială: 1 0 0 cos p cos Ti cos p sin 7 sin p 0 cos £ -sin € 0 \ / cos S cos oc\ ,n £ cos fc sin = 5h43m01s,5 ; p = 45°07'46" ; £ =23°27'26" К e z o 1v a r e : Se aplică formulele lui Gauss triunghiului ■feric -TTP6’ (fig 2 19) (unde a = 90°-S , b = 90°-p) Se obține: sin S = cos £ sin p + sin £ cos p sin'X , cos £ sin oc = - sin £ sin p + cos £ cos p sin 'Л , (2 11) cos £■ cos oc = cos *X cos ți , unu, în forma matricială: / cos § сов ос \ /1 0 0 \ / cos p cos Л\ cos S sin oc j = 0 cos £ -sin & I ! cos p sin 7 j (2 12) \ sin S ) \ 0 sin £ ' cos 6/ \ sin в / Aceste relații permit trecerea de la coordonatele ecliptice la coordonatele ecuatoriale Pentru aplicația numerică, p • (2 16) Soluția -E corespunde răsăritului, iar soluția +H corespunde spusului astrului Timpul sideral al răsăritului și apusului este dat respectiv de: = ос - H , 0a = ос + H Pentru а calcula locul răsăritului și apusului astrului, ее aplică formula cosinusului pentru latura PS = 90°- S a triunghiului sferic PZS: sin S - cos z sin - sin z сое 4* cos A (2 17) Cînd z = 90°, din (2 17) rezultă: cosA = - sin S /сов , (2 1S) relație care este satisfăcută pentru valorile A și ЗбС°-A, ce corespund apusului și, respectiv, răsăritului astrului 5 Formulele (2 16) și (2 18) se mai pot scrie și în forma: cos В = - tg 8 /tg(90°- 0) , >• 'uită că pentru cei doi observatori, 0^ și C2, care obeer- - 52 vă simultan steaua 6 obținem:: z^ ® 4*2 “ & Ș-t * Eăcînd diferența, rezultă formula: Z1 - z2 = ^1 - ^2 • ■ ' Observație : Dacă astrul 6 este chiar Soarele și dacă observația se efectuează în ziua solstițiului de vară, la a-miază, în 0? (Aesuan, situat pe tropicul Racului} 2^=0°) și, în (Alexandria; z^»?0,2) , obținem - ca și Eratostene -n° = = 7°, 2 3 2 Considerînd'Pămîntul drept elipsoid de rotație, ~~~ avînd turtirea £ sa se arate că între această mărime și excentricitatea unei elipse meridiane, e, e~ xietă relația: 5 9 (1 - Ё ) = 1 - e Rezolvare: Turtirea elipsoidului, 6, se definește prin ?ig 3 2 Secțiune meridienă a elipsoidului terestru, formula: ё = (a - b)/a , unde a - eemiaxa ecuatorială și b - aemiaxa polară a-le elipsoidului (fig 3 2} I Excentricitatea elipsei meridiane (e) obținute prin] secțiunea elipsoidului cu un plan meridian oarecare (figJ 3 2) este, dată de formula: ,1 - 53 - Avem, succesiv: 6 = 1 - b/a => Ъ/а = 1 - € ; e = \/l - b2/a2 =-> b/a = \/l - e‘; Prin egalarea celor două expresii pentru b/a și ridicare la pătrat, obținem relația căutată De aici, de exemplu: e ж '/1 - (1 - ё )2 Aplicație numerică : ё = 1/298,3 și e = 0,0818 3-3 | Să ae exprime coordonatele geocentrice ecuatoriale rectangulare ale unui punct al elipsei meridiane ■ M(x,y), cu ajutorul latitudinii geocentrice 4” li e zo 1 va re : Fie observatorul Fig 3 3 Pentru exprimarea coordonatelor x,y în funcție de, latitudinea geocentrică ii Internul de coordonate Txyz , ales în punctul К de rază vectoa re geocentrică și latitudine geocentrică Ч3' ; K' punctul de pe cercul cir cumscris elipsei, corespunzător punctului M de pe elipsă și situat pe perpendiculara dusă prin К pe a-xa mare a elipsei avînd pi ciorul în Q; Ф = M’TQ este latitudinea excentrică a lui M (fig 3 3) Coordonatele rectangulare ale punctului К în în planul elipsei meridi- 54 ane, au e iile: x = p ccs în funcție de 4м Д-ceasta se face prin intermediul utilizării unghiului Ф Din ecuațiile parametrice (cunoscute) ale elipsei (pentru care, în cazul nostru, parametrul t este chiar Ф): x = a cos Ф , у = b sin Ф = a Vl-e* 2 sin Ф , unde e este excentricitatea elipsei, rezultă imediat relație g = '\/x2+y2 = aVl ~ e2sin2 Ф • Relația dintre unghiurile Ф și se poate obține în modul următor: sin = y/g = aVb-e2 sin Ф / (а'Д - e^ein^ Ф ) , de unde se deduce imediat: sin^ Ф = sin” Ч7* / (1 - e^cos2 4” ) Așadar, prin înlocuire, rezultă relația căutată pentru raza vectoare geocentrica a punctului M: p - aVl-e2/ \/l - e2cos2 9' , de unde coordonatele rectangulare x,y în funcție de Ч5' : x = a '/1-e2 cos ф'/ Ѵі - e2cos2 Ч3' » у = a е2сов2 4” ■ - 3 4 j Să se exprime aceleași coordonate rectangulare ca și în problema 3 3, în funcție de latitudinea geodezică 2 (= - notația din manual) - 55 |i • ? о 1 ѵаге : Alegem acum sistemul de coordonate Txz în planul elipsei meridiane c punctului Mîal elipsoidului terestru (T)(fig 3 4) Coordenatele punctului M(x,z) trebuie să satisfacă atît ecuația elipsei CB): х2/аг + z^fl-e^} = 1 cît și ecuația tangentei i л ' (T) trecînd prin M: 'ig 3 4 Pentru exprimarea со- ■rdonatelor x,i în funcție de latitudinea geodezica z = - (1/tgВ) x +c , iHinetanta c trebuind aleasă astfel ca elipsa Câ ) și dreapta ('/ ) să aibă un punct comun Condiția impusă este satisfăcută dacă avem: c = a Vi + ctg2B - , ixui ce rezultă din discriminantul ecuației de gradul al doi-і'чі care determină soluție comună a celor două ecuații Coordonatele rectangulare ale punctului M în funcție de li rezultă acum, prin calcule simple, astfel: x = a coa В/ V ?- ecein2B , z = a(l-e2) sin В / \/l ~ e‘;Ein2B ' nervație : Neglijînd diferența dintre latitudinea geo-■1'1 ca (B) și latitudinea astronomică ( ), deci nefăcînd distincție între cele două latitudini (В = Ч’) , formulele obținuți servesc și pentru exprimarea coordonatelor rectangulare ale - 56 - punctului Ke('i) în funcție de latitudinea astronomică Ч’ : x = а сов Ч’ / Vi - e^sin^cp , z = a(l-e£ ) sin У se exprimă prin formula: X tg 4” = (1-e2) tg tp 2) Această formulă se poate obține și altfel, fără a cunoaște soluțiile problemelor 3 3 și 3 4, așa cuf не indică în manual 3 7 Să не arate că între latitudinea excentrică (Ф) și la-~ titudinea geodezică (B) există relația: tg Ф = Vl-e2 tgB , unde e este excentricitatea -elipsei meridiane R e z 0 1 t a r e : Din ecuațiile ane : parametrice ale elipsei meridiJ x = a cos Ф , у = a\/l-e2 sin Ф și din ecuațiile: — 59 - cos 4” , У = 6 sio 4” , rezultă relația: Vl-e” tg Ч1 = tg Folosind rezultatul din problema precedentă, 3,6, avem: \/l-e2 tg Ф = (1-e2) tgB , da unde: tg = \/l~e2 tg® • observație : Nefăcînd distincție între latitudinea geodezică și latitudinea astronomică (B = ^),are loc și relația: tg f « x/ÎZe2 tg Ф 2АІ Să se demonstreze că diferența dintre latitudinea geodezică (B) și latitudinea geocentrică (4”) se rxprimă prin următoarea relație aproximativă: B- 4>' = 103 132",4 e2sin(2B) , ■înde e este excentricitatea elipsei meridiane Care este max(B - 4” ) și pentru care cerc meridian ere loc ? I1 ‘ z o 1 v a r e : Pentru determinarea diferenței В - 4” , vos plroa de la formula trigonometrică (cunoscută): tg(B - cp) = (tg в - tg ч>’ )/(i + tg в tg ц” ) ijl vom folosi rezultatul din problema 3 6: tg 4>' = (l-e2l tg В Rezultă: tg(B-T') = e2tgB/(l + (1-e2) tg2B) = r «in В cos B/(cos2B + (l-e2)sin2B) = e2sin(2B)/(2(l - e2sin2B)) бо bimitîndu-ne la primul termen din dezvoltarea în serie a funcției 1/(1 - e£sin2B) , deja e" fiind mic (e2 = 0,00673, vezi problema 3 2) și ținînd cont de faptul că diferența В- Ч3' este un unghi mic, ajungem la formula: (B-Ч”)" = юз 132", 4 e2sin(2B) Dacă luăm В= V , formula de mai sus devine: (Ф~ц»)" = ЮЗ 132",4 e2ein(24)) Observație : Nu este greu de văzut că diferența В - 4*' se exprimă și în funcție de cu aceeași aproximație (ne-glijînd termenii începînd cu al doilea (de ordinul de mărime e^) în dezvoltarea în serie a funcției 1/(1 - е2сов2 4” )) : (В - 4” )" = 103 132",4 e2sin(2?' ) Observăm, de asemenea’, că max(tp-4”) atinge 11’40" pentru Ф = 45°Об', după prima formulă, și este aproximativ egal cu 11',5 pentru Ч* = 45°, după a doua formulă 3 9 Să se demonstreze că lungimea arcului de meridian ^(CiOg) dintre punctele 0^(4’- ) și ^(M’g), unde reprezintă latitudinea astronomică, se exprimă prin formula : r ^2 €Co7o2) = aU-e2)J4)i,a-e2Bin2 4 4^ C'1 că funcțiile х*(Ф) și у’(Ч’) sînt continue, formula (3 2) devine: „ *(£*) = J^V^TF2'!^ i >l, înlocuind x' și y’ cu expresiile date de (3 3), se u- ,lunge ușor la formula care trebuia demonstrată 1 -10 | Să se determine proiecția unei poziții a satelitu- lui artificial al Pămîntului pe elipsoidul de re- 1 Cerință, astfel ca dreapta ce unește satelitul cu punctul Ș de proiecție să fie perpendiculară pe planul tangent Ic î — 62 - elinsoid în sunetul de proiecție Să se calculeze coordonatele geodezice ale punctului de proiecție în funcție de coordonatele geocentrice rectangulare ale satelitului Bezolvare : Vom considera fig 3 6, în care: S poziția satelitului artificial al Pă mîntului; Fig 3-6 Pentru determinarea coordonatelor geodezice ele punctului subsatelit satelitului cunoscute; lecție a satelitului elipsoidul de ж,у - coordonatele rectangulare ale lui S în planul meridian al M - punctul de pro (T), avînd geodezice (L,B) ce UT' ’mează să se determine; 8 - latitudinea geo centi'ică (=declinația) a satelitului S în primul rină se determină r și £ după formulele evi dente (fig 3 6): r = x/x2 + z2 , 8 = arctg(z/x) Coordonatele В și H se determină prin aproximații succe- sive în modul următor: In prima aproximație se poate considera ') , fiind excentricitatea elipsei meridiane ( face după formula: H = \/r2 - g2sin2(B- 4” ) - g cos(E - 4” ) • In sfîrșit, se determină valoarea aproximativă a unghiului (i - ) - 1",16 sin(44>) , știind că excentricitatea elipsei meridiane a Pămîntului este e = 0,0818 3 13 | Să se arate că raportul dintre raza geocentrică a u-“ nui punct al elipsei meridiane și semiaxa mare a e-lipsoidului terestru, în a doua aproximație (în același' sens ca în problema 3-12), se exprimă prin formula: E/a « 1 - 0,0033245 Bin(2«P) - 0,000285 ein(4* lP) 3 14 Care sînt dovezile observaționale ale mișcării de rotație a Pămîntului în jurul axei sale ? Răspuns : Pentru mișcarea de rotație a Pămîntului în ju- 1 rul axei polilor, de la vest la est, în sens opus mișcării a-J celor ceasornicului, pentru un observator situat în emisfera j nordică, dispunem de următoarele dovezi directe: 1) Fenomenul devierii spre est a corpurilor de la vertica-, la locului, în căderea liberă, observat pentru prima dată în j| - 65 - »nuj 1804 (înălțimea căderii fiind de 85m , iar devierea ob-Hnrvată de 11,5 mm ); modelarea matematică e fenomenului se i i - e în cadrul mecanicii teoretice ?) Experiența cu;pendulul lui Foucault, bazată pe proprie- ti pendulului de a-și păstra planul în care oscilează, a forat efectuată pentru prima oară în anul 1851 (lungimea pendulului fiind de 67 m , iar masa sarcinii ~ de 28 kg); modela-!•’>» matematică a fenomenului se face în cadrul mecanicii teoretice, arătîndu-ee'că măsura devierii planului este 15°einH’ unde T este latitudinea geografică a locului i) Devierea corpurilor în mișcare liberă pe suprafața Pă-m'i ntului, observată la începutul secolului XIX de către Coliliile (fenomenul se dovedește experimental în cadrul fizi-ili); devierea de la direcție mișcării, în emisfera nordic;;, ni face spre dreepta, iar în emisfera sudică - spre stânge, fenomenul fiind modelat în cadrul mecanicii teoretice (mi ș-■nirn în sistemele de referință mobile, neinerțiale) Acest fenomen are consecințe importante asupra multor fenomene ge-"II zice, ca: direcțiile vânturilor alizee (în țările tropi-auli-: nord-estice în emisfera nordică, sud-estice în emiafe-П imdică), direcțiile ciclonilor (cunoscuți din meteorologie • llind formațiuni cu distribuții barice aproape concentri-■• , avînd presiunea minimă în centru), devierea curenților tu "unici etc 4) Existența forței centrifuge pe suprafața terestră și a inii tril Pămîntului de-a lungul axei sale de rotație ' ) Observații optice (vizuale și fotografice) efectuate in iijutorul aparatelor instalate la bordul sateliților arii- ьб - ficiali ai Pămîntului 3 15 | Care sînt dovezile observaționale ale mișcării de ~~~ ~ revoluție a Pămîntului în jurul Soarelui? Răspuns : 1) Fenomenul paralaxelor stelare, modelat ma- tematic prin elipsele de paralaxă anuală (vezi capitolul IV) 2) Fenomenul aberației luminii, modelat matematic prin e-lipeele de aberație anuală a luminii (capitolul IV) 3) Obв erva ț ie : De fapt, așa cum a observat astronomul V Ambarțumian, coincidența perioadelor de revoluție în jurul Pămîntului ale planetelor exterioare Marte, Jupiter și Saturn (pe epicicluriie corespunzătoare), precum și ale planetelor interioare Mercur și Venus (pe deferentele lor) cu perioada de revoluție a Soarelui (coincidență neexplicată în sistemul geocentric al lumii după Ptolemeu; vezi capitolul VI), constituie o dovadă directă a valabilității sistemului heliocentric al lui Copernic'(vezi cspitclul VI), deoarece cauza acestor mișcări aparente ale planetelor, ca și a mișcărilor aparente ele stelelor "fixe" pe elipsele de paralaxă anuală, este aceeași - mișcarea Pămîntului în jurai Soarelui - 67 - CAPITOLUL IV FENOMENE CARE MODIFICĂ -POZIȚIILE AȘTRILOR PE CER t±l Să- se stabilească o formulă pentru refracția astronomică în cazul modelului sferic concentric pentru atmosfera terestră |l r z o 1 v e r e ; în acest model se presupune că atmosfera te-■ «intră este formată din straturi sferice, concentrice și suficient de subțiri pentru a putea fi considerate omogene Fie două stra- ’!t> «1 1 Schema refracției astronomicei turi avînd indicii de refracție n și n’ =n + dn Notăm cu r și r ’ = r •+ dr distanțele lor față de centrul Pămîntului, a cărnii rază este rQ (fig 4 1) ■ In planul în care se situează ra za de 'umi;‘ă consideram sistemul le coordonate polare (r,6) avînd polul în T și drept axă verticala locului TZ Lumina care vine de le hntru pînă la observatorul situat-In O descrie o curbă din •• stratul atmosferic cu indicele de refracție n delimitea- бе - ză arcul' BA, asimilabil cu segmentul de dreaptă BA , a cărui direcție formează unghiul ț cu verticala locului Bacă i și i'=i+di sînt unghiurile de incidență în A și respectiv B, atunci din triunghiul ATB se observă că unghiul de -refracție în В este i- d6 în baza legilor refracției și a teoremei sinusului din trigonometrie, avem: sin i’/sin(i - de) = n/n' , (4 1) sin(i - d6) / sin i = r/r* (4 2) înmulțind termen cu termen relațiile (4 1) și (4 2), se obțint sin i’/sini = nr/(n*r’) sau n'r’sini’ » n r sin i , d& stg i d-Ѳ (4 7) - 69 - înlocuind (4-7) în (4 5), se obține: d© + di = - tg i (dn/n) (4 6) liln triunghiul AMT se observă că £ = © + i și atunci: d£ = de + di , (4 9) "iu (4 8) și- (4 9) rezultă: d^ » ~tgi(dn/n) , (4 1-0) Hhinltă ecuația diferențială a refracției astronomice Ea re-i"1 ' tntă variația distanței zenitale, £, datorită refracției utrâtul, elementar din vecinătatea punctului A Refracția astronomică este variația totală a distanței • •ni tale, 5, de-a lungul razei de lumină, cînd aceasta stră-t*l> toată grosimea atmosferei Ea este dată de integrala «iii'btlinie: R“ Jcdt = - Jc tg i (dn/n) , (4 11) •i'uiiată de-a lungul traiectoriei luminoase între punctul su-pn'iur de intrare în atmosferă și suprafața Pămîntului Pin (4 4) rezultă: sini = (ro/r)(nQ/n) sin z *1 я lunci (4 11) devine: 4 - 'j„ (ncrosinz/'/n2r2sin2z )(dn/n) , uo v ° Cno H ■ sin z (1/V(nr/(noro))2 - șin2z )(dn/n) , (4 12) йнніііл Integrala refracției Ea poate fi dezvoltată în serie * и puterile impare ale lui tg z : 70 - R = ă tg z + В tg^z + C tg5z + , (4 13) coeficienții A, BrC, ' deducîndu-se fie din ocservații, fie teoretic, pe baza modelelor de atmosfera terestră гг '' ' - - - - - 4 2 ț Să se determine efectul pe care îl are refracția ™”~ asupra coordonatelor ecuatoriale ale unui astru Rezolvare : Fie 6(ol,8) poziția adevărată a unui astru Fig 4-2, Pentru determinarea efectului refracției asupra coordonatelor ecuatoriale ale unui astru pe sfera cerească, iar po- 1 ziția aparentă a a-l cestuia (fig 4 2) 1 In această situației avem (vezi figura):| 6'6'] * Z1 " a x = dz = R (4 141 în plus, avem: шб"п = ( т-й triunghiul sferic Рйб'^, •Imează cu unghiul q din G în care unghiul din 6^ se apro-Aplicînd acestui triunghi Por- Hnilele lui Gauss, obținem: cos = sin Ч3 cos 8д + cos 4* cos 8, cos Hj , sin z, cos q = sin coe 8^ - cos 9 sin 8 cos Hj , (4 16) sin z-, sin q = sin H-j cos cos H, / сов M s Cu notațiile introduse, sistemul (4 16) (4 17) ff e***x- , — devine: " 6(,r) O а -Д \Л^АЛ 4 3J într-o zi, la momentul sideral Ѳ = 18b43m468,5 > s-a observat din Cluj-Napoca (Ч> =46°45'34") un obiect cosmic, determinîndu-i-se coordonatele ecuatoriale aparente : ct] = 12h04m348,5 , S1 = + 52°53'O5n Să se determine corecțiile doT, dS datorate refracției astronomice, știind că dz = R = A tg Zj , unde A=58",3 Indicație: Schema de rezolvare este următoarea: s— л 1) calculul lui Hx: Hj = Ѳ - of = 6®39m12°, 0 £ 99°48'* V 2) calculul lui M cu formulele (4 17); - Qf 3) calculul lui q cu formulele'(4 18); 4) calculul lui z^ cu formulele (4 18); 5) calculul lui d eîml lik ralaxele de înălțime ale Lunii pentru observatorii din G, rl 02, atunci din patrulaterul TO-^LCK rezultă:' în acest mod se obține 3Vff = 57'O2", 7 Avînd în vedere fa; tul oM raza ecuatorială a Pămîntului este Rq = 6378,2 km și că ■•'к* H0/Ac , se obține distanța medie a Lunii față de Pământ; Zijț = 384 403 km Й ’ ** * r v s t i o : Pentru determinarea pare laxelor diurne ale iilor corpuri din sistemul solar se poate folosi tot metod;: *•*- I ! • in 1,11 din care se deduce imediat: • p •= px +p2 « (Rj/AgHeinz^ + sin Zg) •’* и i и 1 , distanța geocentrică a astrului ff este dată de: 78 Д = Rj (sin + sin zB)/p (4 27) Mg 4 6 Pentru determinarea distanței unui astru Aplicație : Presupunem că = 30°, z^ = 45°■ p = 70* EX' primăm mai întîi p în radiani: p = 70-(91/180 )/60 = 73Т/10Ѳ0 Aplicînd acum formula (4 27), se obține: Д6 = ((1 + V5)*1080/(2•73V) ) Rj = 59,276 Rj 4 11 Să se stabilească formulele pentru determinarea corecțiilor de paralaxa diurnă în coordonate e-cuatoriale Rezolv ar e : Fie sistemul de coordonate rectangulare Txyi cu originea în T - centrul Pămîntului,'avînd axa Tz îndrept» tă în direcția polului nord, axa Tx în direcția punctului T (punctul vernal), iar axa Ту fiind situată în planul ecuatorului astfel încît sistemul să fie drept Considerăm un puno1 К de p:1 suprafața Pămîntului și fie ИГО meridianul locului Sistemul rectangular de coordo; ’e Mx’y'z’ are axele parale-: ie cu axele sistemului Txy- Г'- ъ un astru в situat la dla-l - 7S - тніп Д față de T și la distanța Д' față de К (fig 4 7) n, 4 7 Pentru determinarea corecțiilor de paralaxa diurnă й> »еі astru are, în raport cu sistemul rectangular și polar de (■ «♦•donate cu centrul în T și M, respectiv, următoarele coor- йнн*lai (x' y' z*) - coordonate rectangulare topocentrice; (Л' ос* , S') “ coordonate polare topocentrice; (x,y,z) - coordonate rectangulare geocentrice; (A ot, S) - coordonate polare geocentrice IM» și M(R,6,4”) coordonatele rectangulare și iiei P tiv polare ale lui M față de sistemele de coordonate cu IHdimn în T (R fiind raza Pămîntului) Putem scrie: TC = TM- + M? , t>au pa componente: X - £ + x’ , у » T) + y* , Z = î; + z' , ceea ce duce imediat la: х’«х-£ , y’«y-T] , z'=z-ț , (4 28 aude: x » А сов ț cos oc , x' = A'cos 8’cos oc', £ = R cos tf’cos 6 , у = A cos £ sin a , у'® Д'соѳ S’ein , r, «R cos cp’eine , z •= A sin 8 ; z' = Д’ ein £ ’ j £ = R ein 4” | Așadar, formulele (4 28) ee pot scrie sub forma: Д'сов S’coe ce' = Д coa S сов oc - R сов Д'сое(£*-£) ■- Д ~ R(®e® ‘■p’-sin S +cos ‘f'cos 5 coa(0-a)J A’ein( £’-£)=- R(sin 4” coe S- - сое 4” sin S cos (©-«■) • 4' j Din (4 30) și (4 31) rezultă: tg(oc'-oc) - R cos ф'еіп(Ѳ~ ct)/(A cos 8 - R сое Ч”сое(е-ос)и tg( S S) = - R(sin »l>rlmă în unități astronomice Atunci rezultă imediatr A = 206 265 / 0,76 = 271 400 UA Deoarece la l = 63 304 UA, rezultă că: A- (206 265/31")-(1/63 304) = 3,26 /0,76 = 4,35 a l I ■* I 1 Să ee determine corecțiile de paralaxă heliocentri- ! fc- ' „ л л ” I ca, atît in coordonate ecuatoriale, cit și in cdor- j Itbttio I" ecliptice В* * " l v a r e • Se consideră un sistem de axe trirectangular ■"-1,1, Tx'y'z1, cu originea în centrul T al Pămîntului, avînd I»* i• ' îndreptată spre punctul vernal, iar planul Tx'y’ со- ■te1 1 1 o planul ecuatorului Cu centrul în Soare, care es- ■ț иіііігм in planul eclipticii, se consideră și un alt sistem ■ HH frti>ntite, Sxyz, avînd axele paralele cu axele sistemului M*f'- ■ Coordonatele geocentrice ecliptice ale Soarelui sînt Иі \ "> unde a este distanța Pămînt-Soare '■ o stea situată la distanța A' față de Pămînt || iieiiințri Д față de Soare (fig 4 8) h'■ 1 ( respectiv; rezultă că valorile trigonometrice el- funcțiilor sinus și cosinus aplicate coordonatelor geo-ItHlrlon oc' și 5' pot fi aproximate cu valorile corespunza-flM* - ind ele ее aplică coordonatelor heliocentrice oc, și S i formulele (4 37) - (4 38) se scriu sub forma: \ , -A:1 Bk' i« - -7Г вес S (sin ос cos An - cos oc ein cos £ ) , , Л , , Ѳ 6 (4 IA) * ' ‘ (coe oc sin S coe Ag + ein oc ein S ein (Ag cos £ - - COS S sin (Xg sin e ) , tl«l i »|o «'/,intă corecțiile de paralaxă anuală în coordonate , sensul pozitiv 6^ fiind direcția est a paralelului și 6r] direcția nord a meridianului ecliptic în acest sistem de axe, coordonatele rectan-j gulere ale unei poziții a parente (geocentrice) sînt următoarele £ = ("V - ) cos a , ’ (4 41M 4 = ₽ 1 - p - 85 - Ținînd seama și de relațiile (4 40), se obține: ' £ = -0Г Bir (-?>- 9v„) , (4 42 q => -0Г sin S cos(9i - Eliminând parametrul 9i-7\@ între cele două ecuații de *• I bus, rezultă: Ș2/Jt2 + ^2/(ar2 Bin2p) » 1, (4 43) «ere reprezintă ecuația unei elipse cu semiaxele 7Г și % sin p «(•mită elipsa de paralaxă anuală a stelei In relațiile (4-42) numai parametrul variază în de-futiul unui an de la 0° la 360° Mișcarea poziției aparente (*• «lipsa de paralaxă anuală se efectuează în sens direct, ca sl mișcarea Pământului pe orbita sa circumsolară, iar poziția • înrrintă a stelei diferă cu 180° de poziția Pămîntului pe orbi Im ea Elipsa de paralaxă anuală este o dovadă a mișcării de »«»oluție a Pământului în jurul Soarelui 4 i' | Să se determine expresiile corecțiilor de aberație anuală in coordonate ecuatoriale Дй t) 1 f я ț i e : Ele (oc , S ) coordonatele ecuatoriale ale pofti lui adevărate a stelei S și (ct *, £ ’) coordonatele ecuato-H«l ', p*) va descrii o curbă închisă în Juni poziției adevărate S" I Dacă -ee consideră ■ (fig 4 11) meridianul (il paralelul ecliptic care trec prin C și sistem» rectangular de ccordcnJ te cu axa în dl recția est a paralelul» ecliptic, iar St) în di / i « 1 - О ,164 u m , - 90 - le,s‘ = 1е,Ші - 0B,003u m România aflîndu-se în cel de-al doilea fus orar, legăț ra dintre timpul său mediu t și timpul universal TU (tim] mediu de la Greenwich) este dată de relația: - 1 TU = tffi - 2h (51 Dacă într-un anumit loc de pe Pămînt momentului mediu| îi corespunde momentul sideral 6 , iar miezului nopții miji cil (medii) îi corespunde momentul sideral Q , atunci intei valului de timp mediu t- 0 îi corespunde intervalul de tifl sideral G - 6Q Rezultă că: t„ - 0 = (Ѳ - Ѳ ) u (5І m о Г Pentru observatorul din Greenwich, relația (5 7) devine: tffiG - ти »л(ѳй-ѳо(Л , (5І de unde: eG = (1/Л)ти + eoG Legătura dintre diferența de timp pentru două localitl corespunzînd aceluiași moment fizic, și diferența longitudi lor lor geografice este dată de: HA- HB = V LB’' care, scrisă pentru timpul sideral, unde localitatea A евИ Cluj-Wapoca, iar localitatea В este Greenwich, devine: ®Cj “ ®G * LCj L G ’ sau: I 6Cj = ®G + LCj • ‘ | Formulele (5 6), (5 9) și (5 12) rezolva problema prJ 91 - и я n trecerii de la timpul mediu la timpul sideral în cazul 'numeric dat, avem: TU « 17h35“Ma>2 - 2h= 15h35m43e,2 (1/^)TU = (l/^)-15h + (1/Д1)-35т + (l/^)-43s,2 = - 15 • (lhO9® 856) + 35 ■ (lmOOs,164) + 43,2 • (Iе,003) = = '15h02m27e,84 + 35mO5S,74 43B,13 ■= 15h38®16B,91 • Din (5 9) rezultă: 0Q « 15h38m166,91 + lh13m32s, 6 = 16h51m49e, 51 , Im pe baza relației (5 12) obținem: 0 = 16h51m49®, 51 + lh34m23e, 46 = 18h26m12s, 97 -GJ - ———,——— - La Cluj-Sapoca s-a determinat timpul sideral, obser- vi ndu-se trecerea la meridian a unei stele avînd as-ктіпіп dreaptă oc = 5hO5m42B,O3 Care a fost momentul de Mmediu, corespunzător observației, știind că timpul si- la ohTU era 6oG - 2h24m30s, 5 , iar longitudinea localii ‘ții Cluj-Napoca este Lg-j = 1^34т23е,46 ? V 1 vere : Se știe că 6= poziție PZ6"i și PZS'g rezultă în acest caz: sin sinSj+co8t₽ cos S2 cos H2 ; (5 15) cos cp (cos Sj cos H-, - cos d? cos H2) = = sin Ч5 (sin &2 - sin &^) (5 16) l'unfnd condiția cos Ч’/0 (adică cp ^90°) obținem, prin 1*1 t ?1) агевіп((с/а) cos p) + kit - p , (5 1 unde p~arctg(b/a) , keffi ; atunci: 6 = (-l)k arcsin((c/a) сое p) + кГС - p + ot-, (5 1 Discuție : Deosebim, după valorile lui S>, cîteva cazuj 1) 4 = 0° (ecuator); corespunde cazului a2+b2=c2 Atuncț avem arcsin((c/a) cos p) = 0 și ecuația devine o ecuație' cu soluția generală: Ѳ = k5t - p + otj , cu [3= arctg(b/a) și кеЯ( în acest caz, ecuația admite l-oo soluții reale 2) 4=90? (pol); corespunde cazului a2+b2 c2 Acum ecu«J ția are 2-® soluții reale Se aplică în acest caz formuluj ~ 97 • 'ilnută în forme sa generală '■ '• | Cunoscînd direcția meridiane! și declinația Soarelui la un moment dat, să ее determine ora locală, utili iikI pentru aceasta umbra lăsată de un stîlp vertical de lungime 8 jf» o l v a r e : Să considerăm fig 5 2 Din triunghiul drept- unghic COB putem scrie: tg 6 = e/oc Se mai observă că: în triunghiul sferic de poziție PZS (S - Soarele) se cunoaște Sg și latura ZS: ZS = = 90° - hg • Aplicî nd în acest triunghi sferic teorema sinusurilor, se obține: nln Hg/sin Zg « sin Ag/sin(90°- Sg) , ein Hg ж (BinAg/cos Sg) coa C im ii unghiul orar al Soarelui, Hg, este cunoscut, atunci in, 1 iu este dat de relația: • t ж H„ + E + 12h , m Ѳ Hîvî» | t Iii IC п-в notat ecuația timpului - 98 - Un acuzat de comiterea unei crime tă, la o oră cunoscută, aduce ca ției sale o fotografie a ea făcută în alt oraș Pe fotografie se vede umbra de telegraf, atît pe teren, cît și pe poate verifica faptul că' fotografia a : și la ora la care в-a săvîrșit crima ? 5 7 1 într-o zi cunoscu-dovadă a nevinovă-fața unei case din lăsată de un stîlp zidul casei Cum se fost făcută în ziua Situația este reprezentată schematic în fi, 5 3, în care AB este stî: Fig 5 3 Pentru determinarea datei și orei după umbră pul, iar ACI) — umbra stîj pului Prin măsurători b ’ | Oă ее reprezinte grafic „ecuația timpului”: E = t - t , m a •> lueind datele numerice din Anuarul Astronomic pentru a- ичI Discuție • * I rnre eînt scările de timp folosite astăzi în astronomie ? |gf I urin : a) Timpul universal (TU) Dacă se cere o pre-■lllr mnl bună decît 0s,01 , trebuie specificată forma în câți itU loloalt (TUO = timpul solar mediu, obținut direct din jWlț’Hlll netronomioef TUI = TUO corectat de efectele mici-ale Pămîntului relativ la axa sa de rotație; în 'Iii;1 > TUI corectat de efectele micilor fluctuații în viteza de rotație a Pămîntului) bl Timpul efemeridelor (TE), timpul dinamic terestru (TOT) ■ *4m|n>i dinamic baricentric (TDB) Scara TE a fost înlocu-|H 9țti ,»') ^raport cu punctul material (>,,■ se numesc elementele vectoriale ale orbitei, componentele lor reprezentând coainușii directori ai axelor de coordonate or bitale ; P O I cos И сов cu - X ПТ XI g-ln fl ВІПСО cos i, - cosii sinЮ- sinii sin i\ - sinii cos со cos i, py S \pz % Sz = sin П cos CU + + cos £1 sin Ш СОВ i sinii sini, -sin 11 sin w + -cos H sini + cos ii cos (£> cos i, cosii sini, cosi / la care se adaugă și formulele pentru componentele vitezei fată de sistemul de referință, fundamental, Sxyz 7 14 j Să se stabilească succesiunea funcțiilor ce reprezintă soluția generală a problemei celor două corpuri, în cazul mișcării hiperbolice • Kezolvare : Cazul mișcării punctului material (P,m) în raport cu punctul material (S,cM ) pe hiperbolă diferă de cazul mișcării pe elipsă numai prin faptul că, în ecuațiile în care apar excentricitatea și semiaxa mare, trebuie considerat e >1 și a Atunci rezulta ecuația: 6-3/3 = (t - T) (7 10) ,j 2) Cazul fa 0 Introducînd variabila auxiliară E cu ajutorul relației I r - a(l - сов E) , rezultă: f ® j0' /a(l - ccs E)/(l + cos E) а яіп E dE = (t - T) - 141 și, după calculul integralei: E - sin E = n(t - 'T) , (7 11) unde n = a~J/ 3) Cazul h > 0 Deoarece în acest, caz variabila uniformizatoare a rela- ției dintre г ți t este imaginară, putem pune, introducînd argumentul real H : r « |a|(chH-l) , ceea ce dă, în mod analog cazului precedent: sh Я ~ H (7 12) und e у -3/2 Vom reprezenta acum grafic funcția r»r(t) în toate ce- Fig Ț 4 Variația razei vectoare în funcție de timp în cazurile mișcării rectilinii le trei cazuri ale mișcării rectilinii (fig 7 4) în cazul h 0, curba reprezentată prin relațiile obținute are un singur punct în care r = 0 Pentru t-r + oo, respectiva curbă se apropie asimptotic de două drepte simetrice față de dreapta t = T Se observă că, rămînînd în domeniul noțiunilor matematice abstracte, mișcarea este reprezentată, cu ajutorul fermul»' lor respective, pentru toate valorile timpului te(- ) , indiferent de prezența ciocnirilor la momentele reale Cu alte cuvinte, la momentul ciocnirii, punctul (P,m) se comportă ca și cum ar ricoșa de la punctul (S,cM), adică ca și cum a-ceste ciocniri n-ar înceta mișcarea Pentru obținerea vitezei relative a punctului material (P,m), se derivează în raport cu timpul formulele obținute pentru raza vectoare, r 7 17 Să se demonstreze că ecuațiile mișcării relative în — ргйЪІеійа celor două corpuri posedă și următoarele integrale prime (integralele lui Laplace): - 143 у гж - гжу = £ , z г* - r*z = fz , unde am notat, pentru prescurtare (introducând, de fapt, o variabilă auxiliară): • • V « Г = r r S XX + yy + zz , iar f , f , £ sînt constantele de integrare (componentele X у 7j * vectorului lui Laplace, f ) Rezolvare : Derivînd expresia lui r* în raport cu timpul: » - - ~ ? -2 -2 Г = XX + yy + ZZ + X + у + z și folosind ecuațiile mișcării relative și integrala energi- ei, obținem: Г* = - (^/r3)(x2 + y2 + z2) + x2 +y2 + z2 = = -ft/r + 2ул/г + h = ^i/r + h Derivînd și această relație în raport cu timpul și ținînd cont de notația din enunțul problemei, se ajunge la următoarea ecuație diferențială: r* = - (^/r2) r = - (^/r3) r r = - (p/r3) r* , ecuație asemănătoare cu ecuația vectorială a mișcării punctului material (P,m) față de punctul material (S,cA), unde = G(t4t + m) și r =\Zx2+ y2+ z2 Păcînd combinațiile ir*-r*x etc,, obținem următoarele ecuații (asemănătoare întrucâtva cu cele din care au rezultat integralele ariilor): x?- r*x = d(xr*- r*x)/dt = 0 , у rK - rxy = d(y r* - r*y )/dt = 0, z гл - rMz = d(zr*- r*z)/dt = 0, - 144 - de unde, prin integrare, ajungea la integralele lui Laplace din enunțul problemei în problema următoare vom arăta că integralele lui Laplace pot fi puee și sub altă formă, dacă ne folosim de integralele ariilor 7 16 i Să ее arate că integralele lui Laplace pot fi puse sub următoarea formă vectorială: «ф Д,, (2*/r) r + c xr = - f , echivalentă cu următoarele integrale scalare: - (^i/r) x - (cyz - czy) = fx , - (^/r)y - (сгх - cxz) = f , - (^/r) z - (cxy - cyx) = fz •' Ее zo 1 var e : Pornind de la expresia derivatei r* din pro , , p • 2 ' • 2 blema precedentă: r = -д/г + х + y +z și utilizînd și ex-/ preeia lui r* din aceeași problemă, integralele lui Laplace, după calcule simple, se pot scrie sub forma: - x + y(xy - z( zx - xi) * fx , - (р/т} у + z(yz - zy) - x(xy - yx) = fy , - (p/т) z + x(zx — xz) — y(yz - zy) = fz Polosindu-ne de integralele ariilor ale ecuațiilor miș- cării relative: yz - zy ж cx , zi - xz ж cy , (7 13) xy - yx = cz , unde cx, Cy, c2 sînt constantele ariilor, reprezentînd compo- - 145 - rentele vectorului c (vectorul moment cinetic), rezultă imediat integralele scalare din enunțul problemei Heamintindu-ne că integralele ariilor folosite sînt reprezentarea analitică a integralei vectoriale a ariilor: rx?» ?, (7 14) rezultă și forma vectorială a integralelor lui Laplace din e-nunțul problemei 7 19 Să se demonstreze că între constantele de integrare f (f , f,r, f ), F(с , c , c ) și h există relațiile: X J' Z Л У Z f ■? = o , fC = + hCL ; unde f= I f I, c = icl, ceea ce arată că din cele șapte integrale prime ale ecuațiilor mișcării relative, din cadrul problemei celor două corpuri, numai cinci integrale prime eînt independente Не zolvare : Dacă înmulțim scalar cu c ecuația vectorie- 1ă a lui Laplace: cxr+(^/r)r = -f, (7 15) m reamintim faptul că: F • F = C (7 1b) (ecuația planului mișcării) și ținem seama de proprietatea- pro- dusului mixt: c • (c xr) = r • (c xc) = 0 , răzuită: F • F = 0 , (7 17; ndică vectorul lui Laplace se află în planul orbitei mișca: - 146 - relative (ca și vectorul r ) A doua relație se obține prin calculul direct al expresl 2 2 «2 ilor lui f și c , din care se ve elimina r1 : f2 = f2 + f2 + f2 = г2 г*-2+ r*2 (x2 + y2 + z2) - 2г*гж(хх + yy + + zz) = r2( u/r + h)2+гж2(22і/г+h) - 2r*2(^i/r+h) » = + h)2 - r*2 h ; c2 = c2+c2+c| = (yz - zy)^ + (zx - xz)2 + (xy - yx)2 = = (x2+ y2+z2)(x2+y2 + z^) - (xx+yy + zz)2 = = r2 (2?i/r + h) - r*2 Din aceste relații se obține imediat: f2-hc2 = yp2, (7 1K) ceea ce trebuia demonstrat 7 20 Să se demonstreze că orbita relativă a celor două I puncte materiale, (P,m) și (S,c4t), este dată de so- I luția (comună a) următorului sistem de ecuații: CXX 4 Cyy + C2Z x 0 , ^r + fxx + fyy + fzz - o2 = 0 , prima dintre aceste ecuații reprezentînd ecuația planului I orbitei, iar cea de-a doua - ecuația unei suprafețe de or- I dinul al doilea 1 Indicație : Trebuie arătat că punctul mobil P(x, y, z) II află întotdeauna pe suprafața de ordinul al doilea- din enun» țul problemei (în planul orbitei acest punct se află în mol evident) Se pornește de la ecuația vectorială a lui Laplaol satisfăcută de coordonatele rectangulare ale punctului mobil - 147 - 7 21 Să ее demonstreze că direcție semiaxei mari a secțiunii conice din probleme precedentă este arătată de vectorul lui Laplace Indicație : Va trebui demonstrată conjectura că suprafața dată de ecuația yir + fxx + fyy + fzz - c =0 este o suprafață de rotație, a cărei axă de rotație este o dreaptă ce trece prin, originea coordonatelor (centrțșS atractiv, S), vectorul ei de orientare fiind vectorul f I 7 22 Fie un sistem de coordonate rectangular drept, avînd ca plan S^rj - planul orbitei relative, axe S£ fiind orientată după vectorul lui Laplace, f , iar axa -după vectorul moment cinetic c (coordonate orbitale) Să se stabilească relațiile de transformare dintre sistemele de 1 coordonate rectangulare Sțr]£ și Sxyz Rezolvare 7 5 Coordonatele orbitale sînt definite în fig Fig 7 5 Pentru definirea coordonatelor orbitale - 148 - Din faptul ей elf rezultă că vectorul director al axei St) este vectorul c*f Cosinușii directori ei axelor de coor* " donate orbitale SJș , S«j , S£ fa\ă de axele sistemului Sxyz vor fi dați de următorul tabel: n X cx/c У Iy/f ^zfx-cxfz>/ Cy/C z fz/f (cxfy - CyfxV(cf ) cz/c Atunci relațiile de transformare se scriu astfel: X = (fx/f)e + ((oyfz-czfy)/(cf))r1 + (Cx/c)ț , У = (fy/f)’E + ((c2fx - cxfz)/(cf))n + (cy/c)ț , z = (f2/f )£- + ((cxfy-cyfx)/(cf))r] + (cz /c)ț 7 23 Să se exprime componentele vectorului lui Laplace» fx»fy, fz, în funcție de elementele unghiulare ale orbitei fl - longitudinea nodului ascendent, i - înclinarea orbitei și JL «a Z z = (fz7f)rcos v 4- {(cxfy - cyfx)/(cf)) r sin V - 149 - Comparînd aceste relații cu relațiile cunoscute din soluția generală a problemei celor două corpuri (care pot fi citite șl din triunghiurile sferice xPN, yPS și zHf din figura 7 5): x/r >» coeQ ccs(co + v) sinii sin(co + v) cos i = = (cos b> cos Л - ein to sinii cos i) cos v - - (sin cu сое £2 + сое ш sin£1 соя i) sin v , y/r = cos(90°—11) cos(w + v) + sin(90°- П) sin(w + v) cos i » = (cosa) sin H + sin cocos Л cos i) cos v - - (sin a) sin£2 - cos io coeâ coe i) șin v , z/r » coe 90° coe(w+▼) + sin 90° eln((0 + v) cos(90°-i) = » sin as ein i cos v + cos w sin i sin v , rezultă relațiile căutate: fx/f = cos (о сое П - sin to sin £2 coe i , fy/f « сое O) ein £1 + sin o) сое 12 coe i , (7 19) fz/f “ ein Ш sin i , precum și relațiile: (cyfz-czfy)/(cf) ■ -einwcoell - cos o> sinii coe i, (c^f^ - cxfz)/(cf) « - sin W sinii +cos tO/CosQ cos i, (7 20) (cxfy - cyfz)/(cf) ж совы sini 7 24 Să se stabilească ecuația orbitei pe calea rezolvării sistemului de ecuății din problema 7 20 I 12/ѴХг t-y*1 +Z£ + C X - cxy , - 152 - după care, cunoscîndu-ee f = vf| + f+ f| , Q și 1, unghiul se obține din relațiile (7 19) (problema 7 23): f /f = cos O) cos XI - sin 60 ein П cos i , X f /f = cos &j sin O + sin W cos XI cos i , f ,,/î = sin ș e ,r v a ț j,,$ : Trobleaa este rezolvată pentru cazul mișcării pe elipsă ftmt ru eszwrile mișcării pe hiperbolă, parabolă șau în linie dreaptă, feriBu-lele care au sînt comune (valabile pentru orice secțiune conică) ee înlocuiesc cu formulele corespunzătoare de la problemele 7 14, 7 15, 7 16 7 26 Să se arate că: f 20Г - (1/(2%1)\C rdv = b = aVl-e2 , unde r este raza vectoare, iar v - anomalia adevărată sile punctului material ce se mișcă pe o orbită eliptică, de semiaxe a și b și de excentricitate e , în condițiile problemei celor două corpuri Dar din c = Vp®(l-e2) și n Rezolvare : Exprimăm r și dv în funcție de anomalia excentrică, E : r - a(l - e cos E) ; E - e sin E = n(t-T) , de unde rezultă: dt = (1/n) (1 - e сов E) dE , rdv e (e/r) dt = c dt/(a(l - e cos E)) = (с/(na')) dE \/Ді/а3 se deduce imediat: ©/(na) « a\/l-e2 - b în final rezultă: f 2SC r 23t f 2 л r (1/(211)) V rdv = (l/(2Jt))\r (c/(»a))dE= (1/23T))\ aVl-e-Qr J U *3 W - v « a\/l~e^ c b - 154 - Observație: Folosind formulele (6 5) și (6 6): r = -a(l — e cos E;, r cos v = a(cos E — e), obținem nentru cos v expresia: cos v = (cos E — e)/(l— e cos- E), de unde, prin diferențiere, rezultă: dv = (VI — e^/Cl — e cos E)î dE Deci: r dv = a Vl ~ e2 dE 7 27 I Să se arate că în cazul mișcării eliptice: —J rT 2 (l/T)}0 r dt = a(l + eV2) , r fiind raza vectoare, iar T - perioada de revoluție pe elipsă Indicație : A se vedea rezolvarea problemei precedenteJ 7 28 Să se demonstreze că în cazul mișcării eliptice: unde L este lungimea elipsei, iar d& - elementul de arc al elipsei Re z o 1 v a r e : Avem: , C L Г L (l/b)jorâî = (1/L)2j0 a(l - e cos E) d£ = Von: calcula integrala: I = \сов E dt a - (ae/L)\ cos E dfc î o - 155 - Vom folosi coordonatele parametrice ale elipsei: x ~ & cos E , у = b sin E ; dx = -asinEdE, dy = b cos E dE Rezultă, conform formulei lui Cauchy: dt = VTdxJ^+Tdy)2 = \/a2sin2E + b2cos2E dE Atunci integrala I se exprimă astfel: 1 - \QcosE dC = \0 Va2sin2E + b2cos2E cos E dE •+ о2Я + Jot + b2cos2E сов E dE Punînd E = E^ + 3t , a doua integrală devine: C 2JT f r- 5^, va2sin2E + Ь2совгЕ сов E dE = = ~^0\/a2sin2E^ + b2cos2E-^ cos E-^ dE-, , în urma căreia: CL 1 = 30 cos E d£ = 0 și obținem în final: (1/L)^O rdC «■ a Observație : Comparînd valorile medii (integrale) de la problemele 7 26, 7 27 și 7 28, se vede că cea din urmă este a-cea „distanță medie" care figurează în legea a treia a lui Ke-pler 7 29 Să se arate că, în cașul orbitei eliptice, anomalia excentrică, E, se exprimă, în funcție de anomalia mijlocie, M, cu precizia e2, prin formula aproximativă: tg E = sin M /(cos M - e) - 156 - I aii ca ți e : Se folosește ecuația lui Kepler și se cezvol» tă în serii de puteri funcțiile sin M și cos И 7 30 Să ее demonstreze că componentele vitezei punctului material (?,m) după axele' orbitei sale eliptice în raport cu centrul atractiv (S,cM), în cadrul problemei celor două corpuri, sînt: dx/dt = - (c/p) sin v , dy/dt » (c/p)(s +cosv) , unde c este constanta ariilor, p- parametrul orbitei, e-excentricitatea, iar v» anomalia adevărată Să ae scrie formula vitezei I V 1 pentru ve(0, 25V) și pentru vîrfurile elipsei (perimetru, apoastru) Rezol va y e : în sistemul de coordonate Oxy (fig 7 6), 09' ordonatele punctului Г Pig 7 6 Componentele vitezei în mișcarea eliptică donatelor x și у , obține® : dx/dt « (dr/dt) coe v tyfâ t « (dr/di) sin v sînt: x s ae + r cos v , у = r sin v Ținînd seama de faptul că: r = r(t) , V = v(t) , și derivînd în raport ■ timpul expresiile coor ~ r (dv/dt) sin v , + r (dv/dt) cos v Folosind e,cuația orbit si, expresia vitezei radiale, dr/dt - 157 - și integralele ariilor (din problema 7 13, de exemplu), rezultă: dx/dt = (c/p) -e sin v eosy - (e/p)(lcos v) sin v = = (c/p) (e sin v cos v - (1 4 e cos v) sin v) = = - (c/p) sin v , ' (7 24) dy/dt - (c/p)(e sin" r + (1 + e cos v) cos v) » = (c/p)(e +cosv) Viteza V ее obține astfel: din relația V -"Ух2 + У2 , ți- nînd seama de expresiile (7 24), se deduce: V - (c/p)y 1 4 e2 + 2e cos v (7 25) în vîrfurile elipsei (TT, A): ѴП “ V|v=0° “ (c/pK1 + e) ’ vA = vUw = (0/p)cl-e)- 7 31 I Să se determine masa planetei Jupiter cu ajutorul mișcării satelitului său Io, care se mișcă în jurul planetei pe orbită circulară cu perioada T = lz,769 , la distanța a = 421,6 • IO"1 tan Rezolvare : Ne folosim de legea a treia generalizată a lui Kepler, sub forma: T2(«M + m) = 132,7 ’ IO-16 a3, în care perioada T se exprimă în zile mijlocii, distanța a în km, iar masa în mase terestre (a - masa satelitului se neglijează în raport cu masa planetei) Cu valorile date, obținem: - 158 Л =■- 132,7 • ІО-16 а i Soare este o elipsă foarte■ a> agită - 159 - Viteza V, într-un punct al orbitei de rază vectoare r, determinîndu-se pe baza integralei energiei, după formula: V = Va \?2a/r - 1 , unde V este viteza medie corespunzătoare distanței a , semi-axa mare a orbitei eliptice se obține astfel: Vq = VaV2a/q-l (formulă rezultată din cea anterioară pentru r = q); folosind formula pentru viteza circulară pentru cometă (r = a): Va = 29,78/Va (km/s) , ne obține: V = 29,78 V2/q - 1/a , q de unde: a = (29,78)2/((29,78)2• 2/q - V^) = = 886,85/(886,85 • 2/0,587 - (54,52)2) = 18,OUA Excentricitatea orbitei eliptice se calculează conform formulei: q = a (1 - e ) , obținîndu-ee: e = 1 - q/a = 1 - 0,587/18,0 = 0,967 Conform legii a treia (simple) a lui Kepler, obținem pe- ii rada de revoluție T (în ani siderali) după formula: T = aVa , •înde a este semiaxa mare exprimată în UA Se obține: T = 18VÎ8 = 76 ani N "J' ă : Hi*1 calculele de efemeride a rezultat momentul trece-■11 la periheliu, t=9 februarie 1986, la distanța heliocen- - 160 - trică q = 0,567 UA, generația actuală de oameni fiind martoră a apariției pe cer a celebrei comete Halley La Observatorul Astronomic din Cluj-Kapoca, cometa a fost pentru prima datli observată (în prima perioadă de vizibilitate) în seara zilei de 14 noiembrie 1985, în vecinătatea stelei 37 Tauri, prezen-tind o strălucire aparenta de 7,5 - —— -jl 7 33 i Sa se calculeze pentru cometa Ikeya-Seki, avînd da- I tele q = 0,0083 UA, Vq= 480 km/s, aceleași mărimi ca 1 și în problema precedentă • ■ — — —I R ă e p u n s : ^сігс = 327 km/e , Vpar = 461 km/s , adică V este mai mare decît valoarea vitezei parabolice: cometa, a trecut prin periheliu pe orbită hiperbolică și, deci, ea nu se va mai reîntoarce în vecinătatea Soarelui 7 34 Pentru asteroidul Icarus, semiaxa mare a orbitei es^fl te de 1,078 UA, iar excentricitatea ei e = 0,826 Să I se determine: viteza în periheliu, viteza în afeliu, viteza I medie, viteza în punctul avînd anomalia adevărată v=90°, | precum și viteza circulară și viteza parabolică la distan- I țele respective de Soare Re zolvare : Din ecuația orbitei, pentru 0k ek 1 și v • = 90", obținem: r > 1,078 (1 - (0,826)2)/(l+ 0,826 cos 90°) = 0,342 UA ; iar pentru v = 0° (periheliu) și, respectiv, v = 180° (afeliu) găsim următoarele valori: - 161 - Г)ѵ=0о = q = 1,078 (1 -0,826) = 0,188UA, гІѵ=180° = Q= 1,078(1 +0,826) = 1,96BUA Viteza medie a micii planete ее obține după formula: Va = 29,8/1/1,078 = 28,7 km/s , viteza în periheliu, după formula: = 28,7\/l,968/0,188 = 93,0 km/в , viteza în afeliu, după formula: VQ = 28,7/У1,968/0,188 = 8,86 km/s , viteza în punctul avînd anomalia adevărată v = 90°, după formula : Vr ж 28,7 ^2 • 1,078/0,342 - І = 66,1 km/s , viteza circulară la distanța q» 0,188 UA, după formula: V = 29,8/^0,188 = 68,7 km/s , circ,q Inr viteza parabolică la distanța q= 0,188 UA, după formula: V „ = 1,41 • 68,7 = 96,9 km/s par,q Rezultă că V „ 4 Ѵл 4 n , adică asteroidul Ica- С1ГС țQ Ц рИ Х s Ц rus se mișcă, într-adevăr, pe o orbită eliptică De asemenea, avem: V u excentricitatea de 0,250 și perioada de 2^'4 Gm, știind că masa planetei Marte este de 0,107 mase terestre și raza sa - de 3400 km ? Rezolvare : Pentru determinarea semiaxei mari a orbitei satelitului artificial, se aplică următoarea formulă, echivalentă cu formula (7 26) din problema precedentă: a = 331,2 'Улт2 ,, (7 27) obținîndu-ѳе: a = 331,2 0,107 • (16О)2 > 4650 km Distanța apocentrică (apoariu) ee obține astfel: Q = a(l+e) - 4650 • 1,250 = 5810 km , după care rezultă distanța maximă de la suprafața planetei: HQ = Q-R »= 5810- 3400 = 2410 km Viteza circulară a satelitului este: 7 = 631,3\ÂÂ7Z - 631,3 \/o,107/4650 « 3,04 km/s 8 și atunci viteza căutată în apoariu se obține, ca în problema precedentă, astfel: Vq = Va\/(i-e)/(l+e) , de unde rezultă imediat: VQ = 3,04Vo,750/1,250 = 2,35 km/s - 165 - 7 37 Să se determine semiaxele mari ale orbitelor și vi- i tezele medii ale sateliților Io și Callisto ai ple-J netei Jupiter, știind că acești sateliți se mișcă față de j planetă cu perioadele de lz,769, respectiv 16z,689 și masa ! planetei Jupiter este de 318 ori mai mare decît masa Pa- ' mîntului - — — -——— -;—-—; ! Răspuns : 421,7 • 10'" km și 17,3 km/s; 1883 • 10^ km și 8,18 km/s 7 38 [ Să se determine perioadele siderale și vitezele me-dii ale sateliților Mimas și Phoebe ai planetei Saturn, care se mișcă în jurul planetei la distanțele medii 3 3 de 185,4 • 10 km și 12 960 • 10 km, masa planetei Saturn fiind de 95,2 ori mai mare decît masa Pămîntului R ă s p uns : 22^,6 și 14,3 km/s; 55Oz și 1,71 km/s 7 39 Să se determine masa planetei Marte în mase terestre după mișcarea satelitului său Deimos, aflat la distanța medie de 23,5 • IO'* km de planetă, perioada de revoluție în jurul planetei fiind de lz,26 Se mai dau următoarele date spre folosire: perioada Lunii în jurul Pămîntului este de 27^,32 , iar semiaxa mare a orbitei lunare - de 384,4 • 10^ km R а вр uns : 0,107 7 40 Să se determine perioada de revoluție ipotetică a ‘ Lunii presupunînd că masa Pămîntului ar crește de ! 166 4 ori, iar Luna s-ar afla le o distanță medie de 2 ori mai mare decît cea actuală Se știe că perioada actuală de revoluție a Lunii este de 27z,32 Răspuns : 38z, 63 7 41 Cît ar trebui să fie masa Soarelui pentru ca Pămîn-tul să evolueze în jurul lui cu perioada actuală, dar la o distanță de două ori mai mare decît cea actuală ? Cum s-ar modifica în această ipoteză, perioadele de revoluție ale planetelor Marte și Saturn, dacă distanțele lor ar rămîne neschimbate, perioadele actuale de revoluție ale acestor planete fiind de 1,881 ani și de 29,46 ani ? Răspuns : De 8 ori mai mare; 243z și 10,42 ani 7 42 Să se calculeze prima și a doua viteză cosmică pen-“ tru suprafața planetei Jupiter, ca și viteza satelitului său Io, ce se mișcă pe o orbită circulară de rază egală cu 421,6 • IO''km, știind că masa lui Jupiter («A) este de 318 mase terestre, iar raza sa medie (R) - de 10,9 raze terestre Inaicație ; Dacă distanța r se dă în km, iar masa cor* pului central, Л, în mase terestre, atunci considerăm pentru viteza circulară formula: Vr,țrc “ 631 3 VcAt/r (кш/s); dar dn că masa s/t se exprimă în mase terestre, iar distanța r estt dată în raze terestre, atunci ^c-, ro = 7,91 VJM,/r (km/s) = Vî2>- 4 km/s; V = 17,3 km/e - 167 - 7 43 Să ее calculeze viteza circulară și viteza parabolică față de Soare la distanțele medii ale planetelor Venus (a = 0,723 UA), Pămînt (a = 1,00 UA), Jupiter (a = 5,20 UA) și Pluto (a = 39,5 UA) Explicați rezultatele obținute К £ a p u n в : 35,02 și 49,52; 29,78 și 42,11-, 13,06 și 18,47; 4,74 și 6,70 (toate vitezele sînt exprimate în km/e) 7-44 Semiaxa mare și excentricitatea planetei Mercur sînt 0,387 UA și, respectiv, 0,206; aceleași elemente orbitale ale planetei Marte sînt 1,524 UA și 0,093 Să se determine viteza medie, viteza în periheliu și viteza în afeliu ale acestor planete Răspuns: 47,9 , 59,0 și 38,9 ; 24,1 , 26,5 și 21,9 (toa- te vitezele sînt exprimate în km/s) Pe baza parametrilor orbitali ai Pămîntului, să se calculeze masa Soarelui în mase terestre Răspuns : 333 000 7 «46 Satelitul artificial al Lunii „Luna-10" a avut următorii parametri orbitali; = 350 km (înălțimea minimă deasupra suprafeței Lunii), Hq =1017 km (înălțimea maximă deasupra suprafeței Lunii) Se cer: 1) Semiaxa mare a orbitei satelitului (a) 2) Excentricitatea orbitei satelitului (e) - 168 - 3') Semiaxa mică a orbitei satelitului (b) , 4) Parametrul orbitei satelitului y) 5) Demonstrarea formulei Je leagă anomalia adevărată (v) de anomalia excentrică (,S) pentru o poziție oarecare a satelitului: tg(v/2) = ИВо+1)/(Н *I)tg(I/2) , ч Ч 'I unde R este raza medie a Lunii 6) Perioada de revoluție a satelitului în jurul Lunii (T) 7) Să se arate că satelitul artificial al Lunii are o perioadă 'de revoluție de circa 9 ori mai mare decît perioada similară a unui satelit al Pămîntului care se plasează pe o orbită cu aceeași eemiaxă mare ca și orbita satelitului Lunii Pentru calculele numerice se dau datele următoare : В - 1738 km, ffij/nij- » 1/81,5 , TL= 27s,3 , bl= 384 400 km Indicație : Se utilizează formulele: a = (Hq+HQ)/2 + В ; e (RQ-Ra)/(RQ + Hq + 2R) ; b = ; p = 2(E+HQ)(R + Hq)/(HQ+Hq+2R) ; ‘ - T = Tj-Ja/ajOVai/a^ Vl'+ Шу/ш^ H ă s p u n в : a = 2421,5 km; e= 0,134; b = 2398,4 km; p = = 2375,6 km ; T = 2h54mllB, 8 ' 7 47 Bara cosmică „Luna-2" (căzută pe Lună la 4 septembrie 1959) a avut la distanța de 320 000 km de cen- - 165 - trul Pămîntului viteza de 2,31 km/s Ce viteză avea racheta la distanța de 230 km de la suprafața Pămîntului ? Si dau următoarele valori numerice: H = 6370 km (raza medie a 3 2 Pămîntului),^! = 398 600 km /s (parametrul gravitațional al Pămîntului) Indicație : Se aplică integrala energiei: V2 = 2^j/r + h И й в p un в : 11,11 km/e -— - 7 48 0 navă cosmică este plasată pe o orbită eliptică în jurul Pămîntului, le înălțimea de 230 km deasupra suprafeței Pămîntului, cu viteza inițială VQ - 10,93 km/s La momentul pornirii vectorul viteza este perpendicular pe raza vectoare geocentrică, fiind situat în planul orbitei Lunii In cît timp ajunge racheta la orbita Lunii, de rază egală cu 384 000 km ? Rezolvare : Din enunțul problemei deducem că punctul de pornire al rachetei este tocmai perigeul orbitei sale, a cărui distanță este гп= 6370 + 230 = 6600 km Semiaxa mare a orbitei se obține din integrala energiei pentru mișcarea eliptică: = V2 = /M(2/rn - 1/ai , din care se deduce imediat: a = rn/(2 - r^V2/^) , 3 2 unde ~ 398 600 knr/s și Vo = 10,95 km/s ; valoarea semiaxei mari rezultă a fi a = 450 203 km Excentricitatea se determină după formula: e = 1 - rn/a = 0,9853 • - 170 - Mișcarea medie diurnă rezultă din legea a treia (completă) a lui Kepler: n = V/1/a3 = 2,1 • 10 D s Orbita rachetei intersectează orbita Lunii la distanța r = 384 000 km Anomalia excentrică corespunzătoare punctului de intersecție se determină din ecuația orbitei: r = a(l - e cos E) , de unde rezultă: сов E = (a-r)/(ae) ==> E = 1,425 Atunci din ecuația lui Kepler: E - e sin E = n(t - T) rezultă, timpul de parcurs căutat: t-T = 60 ore 7 49 Punctul material P se mișcă în Jurul centrului a-tractiv S pe o orbită eliptică de excentricitate e = 0,5 La un moment dat anomalia adevărată a lui P este de 90° Care este anomalia excentrică la acel moment ? Care este valoarea anomaliei adevărate atunci cînd E = 90° ? Rezolvare : Se aplică formula: tg(v/2) =V(l*e)/(l-e) tg(E/2) și se obțin unghiurile căutate: E » 60°, respectiv v = 120° 7 50 Unei nave cosmice i se imprimă, la înălțimea de 230 km deasupra suprafeței Pămîntului, a doua viteză cosmică în direcția paralelă cu planul orizontului în cît - 171 timp ajunge racheta la orbita Lunii, considerată cerc cu raza egală cu 384 000 km ? Rezolvare : Deoarece rachetei i s-a imprimat viteza ca-rabolică într-o direcție în care ea se va îndepărta de Pămînt, desfășurarea în timp a mișcării este descrisă de ecuație mișcării pe parabolă: ține а к 260/61 UA Se mai poate obține parametrul oi’bitei, după formula: p = a(l-e2) = 4,224 UA , dar acesta nu este, însă, un element orbital'independent de a sau e Celelalte elemente orbitale nu pot fi determinate, din ■« tivul că razele vectoare К șl sînt situate pe o aceea®) dreaptă 7 55 0 navă cosmică interplanetară se mișcă față de Soare sub influența gravitați; raia a acestuia (fără să fie - 175 - \ Influențată de acțiunile gravitaționale ale planetelor), într-un sistem de coordonate Sxyz, sistem drept, al cărui plan fundamental Sxy este planul eclipticii, axa Sx fiind îndreptată apre steaua Regulus (care poate fi considerată, cu precizia cerută de problemă, în planul eclipticii) La trei momente diferite, pentru pozițiile navei e-au obținut următoarele valori ale coordonatelor sferice ecliptice: t3 = 10X1980, 12h TU: = 2,15 • 108km, ■XȚ = 272°, ₽1 = 0°î t2 = 29X1 1980, 12^0: r2 = 1,86 -10Skm, 9\2 = 289°, ₽£ =0°; t3 = 28X11 1980, 12*^0: r3= 1,47•108km, X3 = 303°, p3 = 0° Să se calculeze elementele orbitale a , e și со Rezolvare : Deoarece la fiecare din cele trei momente valoarea latitudinii ecliptice a poziției navei este zero, nava interplanetară se mișcă în planul eclipticii De aceea, din ecuația orbitei în coordonate polare, avem: rk ~ P/Cl+ecosC^^-^)) > к = 1,2,3 Pe baZa acestor relații pot fi deduse următoarele formule pentru determinarea elementelor orbitale e și to : e = (rk+i" rk)/(rkcos(,xk rk+icos(4+i -w)) ’ k =1’2’ tg 0) - - (Â^cos Л-j +Ă2cos X2 +A,cos Л3)/(А3віп 9^ + + AgBin 9i2 +A^sin Л-,) , unde Л , A2> A-, au respectiv expresiile : — r^ (r2 — r3 ), Aq ® r2 (r3 — ), A = r-j(rx — 2) • Se vede că mai întîi se obține u> j apoi e și p , iar se miaxa mare a rezultă, în final, din formula: a = p/(l-e2) - 176 - Din aceste formule, pe baza datelor numerice ale probii» mei, obținem: a = 124',5 • 10^ km, e = 0,727 , «u = 89°,3 7-56 Astrul 6 se află la distanța R de Pămînt, coordo- natele sale sferice geocentrice ecuatoriale fiind oc,S Fie Xg, Yq, Zg coordonatele rectangulare geocentrice ecuatoriale ale Soarelui (axa IX trece prin punctul vernal (T), iar axa OY este situată în planul ecuatorului cerescț fig 7 8) și x, y, z coordonatele rectangulare heliocentrice ecuatoriale ale astrului Б Să se arate că pentru determinarea mărimilor R, oc, S au loc formulele: R cos ol cos S = x + Xg , R sin o cos S ® у + Yg R sin S z + Z„ ® | Rezolvare : Din fig 7 8, avem relația vectorială: Pig 7 8 Pentru stabilirea relațiilor de bază ale calculului de efemeride R = R0 + r, (7 28) care, ținînd seama de componentele celor trei vectori: I R(X, X, Z) , I тЧх, у, z) , este echivalentă ou următoarele relații scalare: X - X@ + x , Y = Y©+y, (7 29) Z = Zq + z 177 - Pe de altă parte, putem scrie: X = R cos S cos Л , Y = R cos S sin a , Z = R ein S J (7 30 Din formulele (7 29) și (7 30) rezultă imediat relațiile căutate observație : Relațiile de mai sus stau la baza calculului de efemeride; ele se complică în funcție de condițiile di feritelor probleme concrete 7 57 Să se arate că dacă originea sistemului de coordona-' te OXYZ se mută în centrul de masă al sistemului de puncte materiale, G, forma ecuațiilor diferențiale ale problemei celor n corpuri, precum și forma integralelor prime nu se schimbă, dar se modifică constantele de integrare, e-le fiind numite în acest caz constante baricentrice Re zolvare : Notînd cu (x^ , Fig 7 9- Mișcarea sistemului de puncte materiale față de centrul de masă, G z, ) coordonatele baricen- trice ale punctului Io de masă m , pe baza relați ei vectoriale: R^s^+r^, i = O,n-l (7 3 Sistemul diferențial de mai sus posedă următoarele zece integrale prime: с/Иț = a-j , cArj - ay, cM t «= a, , (7 35) сМ^=ар+Ьу, m ; i=o ZZ mp-pi - zp ) = c2 , i=e n-1 y~ niyCZ^ -XiZ1) = Cy , (7 37) i=c n-1 ȘZ mpx^-y xp = cz ;■ 1=0 (3 integrale ale ariilor); * (-U) = h (7 34) - 179 - unde Ч = (1/2) ч-4-) (o integrală a energiei) 1=0 Se observă că, dacă relațiile de transformare (7 32) se rtorivează de două ori succesiv în raport cu timpul și se ține oont de integralele mișcării centrului de masă, G, se obțin e-cuațiile diferențiale în sistemul Gxyz: mi d2xi/dt obținem: n-1 + > mixi = 0 1=1 n-1 + = 0 i=l n-1 A + > i=l mizi - 0 n-1 unde cM = i=o Aceste ecuații permit determinarea vitezei Soarelui (a sistemului solar), -[Xo > XQ > Zo} față de centrul de masă al gru pului de stele ales, G Observat ie : Dacă grupul de stele este grupul stelelor strălucitoare (pînă la magnitudinea■6)) pentru care se dispune astăzi de cele mai complexe date privind masele și vitezele radiale, se obține pentru mișcarea sistemului solar rezultatul că acesta se mișcă în direcția determinată de următoarele coordonate ecuatoriale: « = 270°, € = 4-30°, cu o viteză de 19,5 I km/ s 7*59 j Să se deducă următoarea ecuație diferențială în cadrul problemei celor n corpuri (formula lui La-grange-Je cobi): d2J/dt2 = 2U + 4h , - 181 unde U eete funcția de forță, h - constanta energiei, iar J - momentul de inerție față de originea sistemului de coordonate, origine situată în centrul de masă al sistemului de n puncte materiale II e z o 1 v a r e ; întrucât funcția de forță U, avînd expresia ca în problema precedentă, este o funcție omogenă de gradul (-1) în raport cu coordonatele, adică: U(txQ, ty0, tzo; txx, tylt tz^; ; tx^, ty^, tz^) = “ t-1U(xo,yo, z0; x^y^z ; ; Xn j , У^, ) pentru orice punct al domeniului său de definiție și orice valoare reală a lui t , atunci după teorema lui Euler referitoare la funcțiile omogene are loc egalitatea: n-1 y~~ (Ой/Эх±)х1 4- (эи/3уі)у1+ (au/azi)zi) = (-i)u , i=o în fiecare punct al domeniului de definiție al funcției U, care este, conform cerinței teoremei, și diferențiabilă Folosind ecuațiile mișcării față de centrul de masă, G, membrul stîng al relației noastre devine: У " т^(х^ d2x^/dt2 + y^ d2yi/dt2 + z^ d2z ! /dt2) = - U i=o Ținînd seama de integrala energiei și de expresia energiei cinetice a sistemului de puncte materiale, obținem: У~~ d2x±/dt2 + (dxj/dt)2 + yt d2yt/dt2 + (dy^dt)? + i=o + z± d2zt/dt2 + (dz±/dt)2) = U + 2h Reamintind faptul că momentul de inerție al sistemului de puncte materiale Po(mo), ’ * * * ’ ^n-l^n-l de origi- nea G a sistemului de coordonate Gxyz este: 182 Ц“”Л P P P j = T~" ^xi+ уi+ si > i=o se vede ușor, după două derivări succesive în raport cu timpul ale acestei relații, că membrul întîi al relației de mai ' 2 p sus'este egal cu (l/2)d J/dt" , de unde ee obține ecuația cerută în enunț Observație : Formula lui Lagrange-Jacobi este utilizată în astronomie la cele mai diverse cercetări calitative (studiul stabilității mișcării, al ciocnirilor etc ) I 7 60 Să se deducă integralele prime ale ecuațiilor mișcă-( Г~~ rii relative (față de punctul material Po) ale pro-[ blemei celor n corpuri | Indicație : Se folosesc relațiile de transformare a ecu»« țiilor mișcării relative (față de P ) din ecuațiile mișcăriii absolute, precum și integralele prime ale ecuațiilor mișcării absolute Răspuns : Integralele prime căutate sînt: n-1 n-l # n-1 n-1 - (iam)( OțiA XZ>3-i ■■ ZZmiz-, + 1=1 1=1 1=1 3=X n-1 + ȘZ = cx, i—1 n-1 n-1 # n-1 n-J - (1AA)( J m z - У~m x 2 m z ) + j=l a J 1=1 J 3 1=1 J d 1=1 n-1 +z: e cy > a-l n-1 n-1 n-1 - (1/л)( ȘZ>1X122«1У1" £>1*1 £>?15 + i«i i®i i=i i-i - івз n-1 + EZ = ce ; n-1 , n-1 - n-1 „ - (l/(2cA-)) (( У m^x^) + ( 5 НЦУ*)* + ( S m^z^)'’) + 1=1 1=1 * 1=1 p—, c , p p + (V2) 5 a (x7 +y? + zf) = U + h , 1=1 n-1 unde «Д, = ș» m4 i=o 7 61 j Să se determine condițiile în care mișcările relative, în problema celor trei corpuri, au loc pe orbite I kepleriene (satisfăcînd ecuațiile mișcării din cadrul pro-i blemei a două corpuri)(soluțiile lui Euler și Lagrange) Re z- o 1 t a r e : Pornim de la ecuațiile mișcării absolute din cadrul problemei celor trei corpuri, în care vom folosi indicii 1, 2, 3 (în loc de O, 1, 2 - notațiile din manual) Astfel, în sistemul inerțial OXYZ (fig 7 10) avem ecuațiile: d2R1/dt2 = Ginig/r^j) r^2 + G(m^/r?^)Î4^ , d2R2/dV' = GCm^/rl^) + G(m-,/r| , ) r^ , c~"~'- /df': = G(m,/r^3)F + G(m2/rj2)ij2 • (7 40) Presupunînd că mișcările relative sînt descrise de vecto- rii r12, r2- Și ținînd cont de relațiile: r12 r23 r31 R2 W’- R-, - r2 , -IU , № obținem următorul sistem diferențial: 134 - d^r^/dt^ - -G(mn-нс2) r^/r^g + ^з^з/^З+^31/r31^ ’ р л 2 —* 3 -■> 3 -л- 3 d r^^/dt = -G(nu-Hi^) + Gm^ ^r31/r31 + r12/r12J ’ ^-411 d2?31/dt“ = -G(m3+ml)?31/r31+ Gm2(?12/l12+?23/r23) • Pig 7 10 Configurația generală a celor' trei corpuri P-^, P2, P^ Conform ipotezei, mișcările relative ale punctelor materiale P-j, ?2» au l°c Pe orbite kepleriene, adică sînt descrise de ecuațiile diferențiale-de forma; d2?i2/dt2 = - G % г^2/з^2 , d2?2-j/dt2 = - G %2 r^-j/Tgj > d£rr3;l/dt2 - - G X3 r3i/r3i > (7 42) unde mărimile X-^, X2> X3 reprezintă constante independente de timp 185 Astfel, condițiile ca ecuațiile (7 41) aă fie de for;-,; (7 42) sînt: (m1 +m2 - Xx) r12/r^2 - Из(г23/і23 + ^зАзР = 0 » (ПІ2 + ш^ - %2) ry/r^ ~ ml^r31^r31 * r12/x12^ ° 0 • (m3 + ml ’ X3} 731/r31 “ ffi2 Sâ2Lll £Î5d vecțorii r‘12л г2зх ?31 ®121 ££1^£45Гг Vom considera, de exemplu, că punctele 1, , P^, F- : ■ situate în ordinea din fig 7 11 Avem în acest caz relația e videntă: r31 =’r12 + r23 • ‘7,46 ' Introducînd notația următoare: “ r23/r12 ’ - 186 - obținem relațiile: Г23 = Г12 ’ (7 48) r31 “ + r12 ’ Ținînd cont de (7 48), ca și de faptul că vectorul are sena opua vectorilor §1 г23’ ^in relațiile (7 45) cantitățile X^, %2, X^ pot fi exprimate astfel: %1 = +m2 -т3(1 + 2Л)/(^2(1 + (X)2) , X2 = m2+m3-m1913(2+7i)/(l+ 'A)2 , ' (7 49) Xj = m3 + mj - m2 (1 + A2) (1 + X”2) Deoarece mărimile Xi și m^ (i=1,2,3) sînt constante, rezultă din relațiile (7 49) că și cantitatea > este constantă Aceasta arată că configurația sistemului celor trei punc- te materiale rămîne în timpul mișcării asemenea cu sine însăși Pentru determinarea parametrului ‘Л , să efectuăm înlocuirea re- lațiilor (7 48) în ecu- ațiile (7 42) Atunci obținem: Pig 7 11 Configurația coliniară a celor trei corpuri x2 x3 = X3 , = X1(1+'Ă)3 (7 50) Dacă în una dintre aceste relații înlocuim cantitățile %1,%2,X3 cu expresiile date de (7 49), obținem pentru parametrul SX următoarea ecuație de gradul al cincilea: (mj + m2 )*A3 + (Зпц + 2m2) Л4 + (Зт^ + mp )%3 - (Зт^ + m2) - - 187 - - (2m2 + Зт-j )9i - m2 - m3 = О (7 51) Întrucît în relația de definiție (7 47) cantitatea 'Л ее presupune pozitivă, căutăm numai rădăcinile pozitive ale acestei ecuații Cum coeficienții ecuației (7 51) prezintă o singură variație de semn, rezultă că această ecuație are numai o rădăcină pozitivă (Descartes) Deci punctele materiale P^, P2, P^ se pot așeza într-uh singur mod, în această ordine, pe o dreaptă, pentru ca acestea, pornind cu o viteză inițială corespunzătoare, să descrie orbite kepleriene Mișcările sînt descrise de ecuațiile (7 42) în care valorile coeficienților %2, se determină din relațiile (7 49), după ce se cunoaște valoarea lui Rezultate similare se obțin și atunci cînd cele 3 puncte materiale sînt așezate într-o altă ordine în acest fel, am ajuns la soluția lui Euler a problemei celor trei corpuri 2) Ca2Hl °Î5â Se £2y±ează țoți coeficienții șișțeinului de е£НёІЦ (7д43) Fig 7 12 Configurații echilaterale ale celor trei corpuri Din aceasta condiție rezultă imediat relațiile: r12 = г23 = Г31 ’ X, = x2 = X3 = Aceste relații arată că cele trei puncte materiale sînt așezate întotdeauna în vîrfurile unui triunghi 18В echilateral (fig 7 12) Fixînd într-un plan pozițiile a doul corpuri, cel de-al treilea corp ee poate așeza în două poziții în același plan, pentru a forma cu celelalte două corpuri un triunghi echilateral (fig 7 12) în acest fel am ajuns la soluția lui Lagrange a problemei celor trei corpuri, cînd cele trei corpuri se mișcă în- I tr-un plan fix, pe orbite kepleriene, formînd tot timpul un triunghi echilateral Fig 7 13 Cazul lui Euler cînd cele trei puncte materiale ee mișcă pe elipse Fig 7 14 Cazul lui Lagrange cînd cele trei puncte materiale se mișcă pe elipse - 183 Observație : Soluțiile lui Euler și Lagrange, în cașurile cînd orbitele kepleriene sînt elipse, sînt primele orbite periodice descoperite în cadrul problemei celor trei corpuri (fig 7 13 și 7 14) 7 62 Să se stabilească forma simetrică a ecuațiilor miș- • cării relative în cadrul problemei celor n coi— j puri,prin următoarea alegere a coordonatelor relative (coordonatele lui Jacobi; 1842): - se consideră un sistem de axe paralele cu axele sistemului absolut, avînd centrul în punctul Po: coordonatele relative ale punctului P j(x, y^, z^) se determină față de acest sistem de coordonate;- - prin centrul de masă G^, al punctelor PQ și P^, se duc axe paralele cu axele sistemului precedent: poziția punctului Pg se determină față de acest sistem prin coordonatele relative Xg, уg, Zg; - prin centrul de masă Gg, al punctelor PQ, P, și P9, se | duc iarăși axe paralele cu axele sistemului precedent și i coordonatele relative ale punctului P-j în acest sistem se notează cu x-^, y^ , z^ ș a m d (fig 7 15) Cu alte cuvinte: sistemele de coordonate succesive se aleg astfel ca fiecare punct de masă se raportează față de centrul de inerție al tuturor punctelor prece- dente de mase mo, m mg m , n-1 Rezolvare: de unde: Notăm G^(^i,'i}^, astfel că: » i -> k=;C (7 52) І ‘ ZZ”Vk • etc , k=o unde prin s-a notat suma: - 190 - Potrivit eu definițiile din enunțul problemei: (unde i>l), din care, după înmulțirea cu A și utilizare» relației (7 53), rezultă: 1-1 eALj i 7X^ ~ i—1 ^i-i = ^'i— 7^1 ~ 1 ®tc \7 57) Pig 7 15 Pentru definire® coordonatelor lui Jacobi Pentru a exprima vechile coordonate X, X, Z prin noile coordonate x, y, ®, observăm că, pe baza relației (7 53): ț,- - 5 etC- ’ (7 58) e&u: - 191 ^5 "* 11—1 ж '**^1 ”'£^i —1 ^Z,ț ’ ®tc« (7 59) Prin urmare, pe bana relației (7 57), putem scrie: tAU (X ” Xi+P “Лі-1(Хі - xp ~ (&Ri “ ’ de unde: Zi+l Zi = xi+l eMi-l ^i1 xi ’ etc’ Adunînd membru cu membru aceste egalități, pentru (7 60) valori- le succesive ale indicelui, obținem: Zi " Zo ' X1 ’ Zi+l“Zo =xi+l +mixi/Jti + ті-1хі-1/сЛі-1 + “ + т1х1/еЛ1- (7 61) Derivînd de două ori în raport cu timpul relațiile (7 57) și folosind ecuațiile diferențiale ale mișcării absolute, obținem: «Mi i d2xi/dt2 » (еЛ(і 1/ті) Эи/ЗХ 5 - -эи/ЭХ0 - - Эи/ЭХ - , - Эи/ЭХі 1, etc (7 62) Pe de altă parte, tot din (7 57), avem: Эи/ЭХ0 = - Эи/Эх1 - (m0XA1) эи/Эх2 - - (в0/Л2) Эи/Эх3 - - (т0/еМп 2)эи/Эхп 1 au/ax-j- = эи/Эх-L - (а^/сАЦ) ÎHJ/Sx2 - - (n^AMg) ' 0 După transformări simple, nunț, rezultă ecuația: rază f> = (K/m^/Cnig- m^a 01 care este ecuația unui cerc de centrul 0 situat pe dreapta ce trece prin centrele atractiv» P, și Po, la distanța | x | = (m,/(m?-m, ))a de centrul atrac- 1 g x c x- tiv P- în partea opusă centrului atractiv P2 (care are таен , 1 ’ mai mare) Cu folosirea raportului q = , raza și centrul echi gravisferei se exprimă prin formulele: 6g = - 9) ’ IXgl = aq/(l - q) (7 64) Observație : Demonstrarea teoremei de nsai bus se poate I face și prin următoarea metodă (sintetică) simplă: din condiția pusă rezultă И P,P||/ІІРрР || = \/q ( = conet ) fixe" - problemă integrabilă a mecanicii cerești} Re zoi va re : Presupunem că vectorul viteză inițială a Fig 7 18 Pentru problema elor două centre fixe Coordonate eliptice punctului P se găsește în planul determinat de punctele %’ Pl’ P2 (Unde Po ' poziția inițială a punctului P), саге se consideră drept plan de coordonate,planul Oxy; originea 0 ее ia la mijlocul segmentu- - 19В - lui Р^Т2, ахв Ox “ direcția dreptei ce unește punctele și P2 (de la P2 epre E^) (fig 7 18) Dacă notăm distanța dintre centrele fixe cu 2c , atunci razele vectoare ale punctelor 1-^ și ?2 se determină după formulele : r, = \/(x-c)2 + y2 , - (7 65) r2 =v (x+c)2 + y2 Ecuațiile diferențiale ale mișcării punctului P, atras d« centrele atractive P, și î2 după legea lui newton, se scriu: * x = ЭП/Эх , (7-66) у = эи/Эу , unde funcția de forță a problemei este: U = Gfm-j/r^ +m2/r2) , (7 67) G fiind constanta gravitațională; r, și r2 au expresiile mal sus stabilite Acest sistem diferențial este echivalent cu următorul sistem canonic: dx/dt = ЭЕ/Эх , dx/dt = - ЭН/Эх , (7 6Ѳ) dy/dt = йН/Эу ; dy/dt = - ЭН/Эу , unde: H « T - U , (7 69) T=(x2+y2) fiind energia cinetică Sistemul canonic dispune de integrala primă: H(x, y, x, y) = h , (7 70) h fiind constantă Sistemul (7 66) sau (7 68) neputîndu-se integre direct, - 199 - Euler, Legendre și Lagrange au arătat că este posibilă introducerea unor variabile noi, astfel încît integrala generală a ecuațiilor transformate să se obțină prin cuadraturi Aceste variabile adimensionale sînt determinate prin următoarele formule: %= (гх +r2)/(2c), ул > (гг - r2)/(2c) (7 71) și sînt numite coordonatele eliptice ale punctului din planul Oxy • Este clar că locul geometric al punctelor pentru care “Л= = const este o elipsă avînd focarele în P^ și P2 și semiaxa mare egală cu 2$\c, iar locul geometric al punctelor pentru care yi^congt este o hiperbolă cu aceleași focare și semiaxa reală egală cu 2^ic Sistemul de referință pentru coordonatele eliptice este deci construit dintr-o familie de elipse și hiperbole homocentrice, avînd focarele în centrele atractive P^ și P2 (fig 7 18) Deoarece avem inegalitățile: r-, +r2 > 2c , |rȚ-r2)4 2c, (7 72) coordonatele eliptice satisfac următoarele inegalități: liîii®, -liyiil (7 73) Din formulele (7 65) și (7 72) rezultă: x « -cdiyu , у = c'/t'X2 - 1)(1 -^i2) , (7 74) de unde: x = - c( u'i + 'X P) , * £, f ț У = C(V(1 -/)7^“- 1) *4 - 1 )/(l -yu2)^) Considerînd (X și u drept noile variabile canonice 200 - definind impulsurile generalizate corespunzătoare prin formu lele: к ЭТ/ЭЛ = C2(CX2-^U2)/^2- !))^ , - ЭТ/ЭД = О2(СЛ2 ~2л2)/(1 , (7 76) ecuațiile diferențiale în coordonatele eliptice vor fi d'X/dt = 9Н/ЭЛ' d^’/dt = -ЭН/ЭЛ d^i/dt ЭН/Эр' d^u’/dt = - ЭН/Э^і , (7 77) unde: H = T - U , T U = (с2/2)(^2-^2)(г2/((\2-і)+yl2/(i-^2)) = = (l/(2c2('X2-212)))((ГЛ2 - 1)9\’2 + (1-ул2)^’2) , = (G/(c(9i2-^^((m^+iUg)^ - (m1-m2)yu) (7 70) Sistemul canonic (7 77), satiefăcînd condițiile teoremei lui liouville (precum și ale lui Stăckel), se integrează pria cuadraturi Vom efectua această integrare prin metoda lui Hamilton-Jacobi Ecuația diferențială cu derivate parțiale a lui Hamll» ton-Jacobi este în cazul nostru: (l/(2c2(5\2 -^2)))((-Ă2 - I)(aw/39O2+ (1 -^i2)(-3W/3pt)2) = = (G/(c('A2-^12))) ((m + m2 - (m^-m2+ h (7 79) Căutînd soluția acestei ecuații sub forma:’ W = W1('A)+W2(^i) , - (7 80) ecuația (7 79) va fi verificată dacă punem: O 2 - l)(dW1/d'Â)2 - 2c2h^2 - 2Gc(m1 + m2)'X - 2 2('/M(/i)) 1dg = Vă (t - y) , \^L(Â)/(A2 - 1) = %' , \/2M(g)/^2 - 1) = /U' , unde am notat, pentru prescurtare: ЪСХ) = (V - l)(c2h%2 + Gc(m, + m?)'X + oc) , (7 -86) M(^u) = (^u2 - l)(c2h^i2 + Gc(m1- m2 )yi + oc) Ecuațiile (7 84) permit 'determinarea variabilele1" 7, ț, i p ca funcții de timp, t, și de patru constante arbitrare, 1 oc, după care ecuațiile (7 76) și (7 85) vor da derivate le în raport cu timpul ale acestor variabile: - 202 - i -\/2L(7j/(c2ft2 - u2)) , , 2 2 2 (7>87) jx У2М(^1)/( с2 (7Г -)) în sfîrșit, coordonatele dreptunghiulare x, у și componentele vitezei x, у ale punctului P se obțin după formulele (7 74) și (7 75) Dar acest procedeu de obținere a soluției problemei puse nu este eficient, datorită faptului că Л și se determină dintr-un sistem de două ecuații transcendente cu două necunos cute De aceea, E A Grebenikov, E P Akeionov și V G Diomin I au propus un nou procedeu pentru determinarea lui X și V- , bazat pe folosirea formulelor (7 87) - reprezentînd de fapt integrale prime ale ecuațiilor (7 77) - și pe introducerea unei j noi variabile independente, T, în locul lui t , după formulai \/2dt = c2(^2-^i2)dT (7 8Ѳ) Deoarece, conform formulelor (7 71): c2CX2-^2) = rxr2 , (7 89) noua variabilă crește monoton cu timpul și poate juca ace- lași rol ca și t Ecuațiile (7 87) în noua variabilă sînt: d'X/d’C » ѴЬ(Л) , - (7 90) d^i/dt =Vm(^i) , care ne conduc, după integrare, la următoarele cuadraturi: V (y/bCxi/d-X » T- , w «0 ° ^(vi^))"1^ = T- To , 5- unde (Хо și p Q sînt valorile coordonatelor eliptice cores- - 203 - punzătoare valorii inițiale a variabilei T (=T0), care poate fi considerată corespunzătoare epocii inițiale t ; de aceea Xo și 2% Pot considerate drept valorile inițiale ale variabilelor ’X și ^i Transformarea integralelor eliptice (7 91) - cere depinde, la rîndul său, de rădăcinile polinoamelor (7 86) - ne dă soluția: » Л(Т-т , % h, 06) , (7 92) Л =Л(т"то’ A’ h> 00 ’ iar, după aceasta, formula (7 88) dă și timpul: V2 (t-t0) = c2^ CX2 -Д") dT (7 93) Exprimarea constantelor arbitrare prin valorile inițiale хо’Уо;іо’Уо 0 Pe seama cititorului Răspuns : h » (x2 + y2)/2- а(т3/\/(х^- c)2+y2 + + т2/У(х0+ c)2+y2) ; = Ф(-х0, Д, h, Л) } Ѵг (to-T) » Y(%\ 2*0 h <x) , prin urmare, soluția generală a problemei conține numai patru constante arbitrare independente - 204 Observații : 1) Probleme, celor două centre fixe este re- zolvată și în cazul generalizat în care punctul material P v mișcă în spațiu, iar mărimile c, , c^ și m-j , nip sînt complex» conjugate; în acest caz ее folosesc așa-numitele coordonate «■ lipsoidale 2) Interesul față de această „problemă restrînsă" clasioft a crescut mult în ultimul timp, în legătură cu studiul mișofl rii sateliților artificiali ai Pămîntului și a navelor cosmice interplanetare Astfel, în cadrul acestui model matematic se obține o aproximare suficient de bună a cîmpului gravita^ onal real al Pămîntului, iar ecuațiile mișcării punctului nm terial în cîmpul respectiv se integrează exact, soluția obținută sub formă de cuadraturi eliptice fiind mai aproape de r« alitate decît aproximația kepleriană (în cadrul căreia Pămîn tul este considerat drept punct material) De asemenea, mițo» rea anei nave cosmice ce pleacă, de exemplu, de la Pămînt ețir» Tună, Venue sau Marte poate fi descrisă cu suficientă aproxl maț ie pe baza acestui model, întrucît durata zborului de 1« W corp ceresc la altul este mal scurtă decît perioada de revolu ție a corpurilor respective în jurul „corpurilor lor" centra le 3) Pe baza integralei energiei și a altor formule ale mod| lului matematic al problemei celor două centre fixe - s tît 1h cazul plan, cît și în cazul generalizat (al mișcării în sp«« țiu) - , se pot efectua numeroase cercetări calitative авирги domeniilor și formelor mișcăriilor posibile ale punctului mut» rial P în raport cu sistemul de referință considerat - 205 - 7 69 Două puncte materiale, (К,т,) și (P2,m2) se mișcă într-un plan fix, pe baza legii atracției universale, descriind cercuri în jurul centrului lor comun de masă (G) Un al treilea punct material:, (P, m), „de masă neglijabilă" (în sensul restricției precizate în enunțul problemei 7 66), se mișcă în spațiu, sub acțiunea gravitațională a primelor două puncte materiale, (Грир și (P2,m2) Să se studieze mișcările acestui sistem de trei puncte materiale, atît cantitativ cît și calitativ („problema re-strînsă a celor trei corpuri") Indicații : Se vor scrie, mai întîi, ecuațiile mișcării „maselor finite" m^ și mo în sistemul de coordonate GXYZ (sistem fix), care admit soluția particulară indicată în e-nunțul problemei Se vor scrie, aprsâ» ecuațiile mișcării absolute a „masei neglijabile" în rapsnrrt, c® același sistem fix de coordonate, GXYZ După aceasta se те tfege un sistem de axe de coordonate mobile, Gxyz, ccnusiderâid Gm2 drept axă Gx, iar axa Gy perpendiculară pe direcția m^nu, , efectuîn-du-ѳе apoi transformarea ecuațiilor mișcării din sistemul fix GXYZ în sistemul mobil Gxyz Se va arăta că ecuațiile diferențiale ale mișcării, atît în sistemul de referință fix, cît și în sistemul de referință mobil, admit o integrală primă a-semănătoare integralei energiei (integrala lui Jacobi) In cel de-al doilea sistem de coordonate, pe baza integralei prime obținute, se va studia cazul V=0 (construirea „suprafețelor de viteză relativă zero", nuntite suprafețele lui Hill, pentru valori date ale constantei de integrare, h) și se vor delimita regiunile din spațiu (în secțiunile planelor de coordonate) în care mișcarea punctului material (P, m) este posibi- lă Pentru simplitate, ѳе та alege sistemul de unități de ml-sură astfel ca, pentru constanta atracției universale, să a-vem G = 1, Observație : Cazul problemei restrînse a trei corpuri »• re o bibliografie deosebit de vastă, jucînd un rol important în cercetările asupra soluțiilor periodice ale unor sisteme dinamice, precum și asupra mișcării Lunii, sateliților artificiali, sistemelor stelare binare etc 